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PREFACE 


This  work,  has  two  purposes:  (1)  to  provide  an  introduction  to  the 

fundamentals  of  fracture  mechanics,  and  (2)  to  provide  a convenient 
source  of  fracture  design  information.  The  presentation,  limited  to 
isotropic,  homogenuous  materials,  is  aimed  at  practicing  engineers  rather 
than  those  engaged  in  fracture  mechanics  research.  In  the  past  20  years, 
a great  volume  of  material  has  been  written  on  research  results  in 
fracture  mechanics  but  it  seems  that  there  has  been  too  little  effort 
devoted  to  distilling  and  summarizing  this  information  and  in  dissemi- 
nating it  to  designers.  Consequently,  engineers,  attempting  to  gain  an 
overview  of  the  subject,  are  frequently  lead  into  a bewildering  maze 
of  technical  papers,  ASTM  special  technical  publications,  proceedings 
of  conferences,  and  the  like.  Practicing  designers  have  little  time 
for  this;  therefore,  an  effort  has  been  made  here  to  present  the  funda- 
mentals of  fracture  mechanics  in  compact  form  in  a single  volume.  At 
the  same  time  a substantial,  but  edited,  compendium  of  solutions  for 
fracture  analysis  has  been  included. 

Chapters  1,  2,  and  3 cover  background  information  dealing  with  the 
history,  philosophy,  fundamental  concepts,  and  basic  theory  needed  for 
a general  understanding  of  fracture  mechanics.  The  effects  of  the 
plastic  zone  in  the  region  of  the  crack  tip  are  considered  in  Chapter  4. 
Techniques  for  both  plane  stress  and  plane  strain  fracture  toughness 
testing  are  briefly  discussed  in  Chapter  5.  Chapter  6 contains  several 
example  solutions  to  some  commonly  encountered  problems  in  the  applica- 
tion of  fracture  mechanics.  Solutions  for  the  stress  intensity  factors 
for  several  different  types  of  problems  are  compiled  in  Chapter  7.  The 
stress  analyst  or  engineer  who  already  has  a fundamental  knowledge  of 
fracture  mechanics  should  find  Chapters  6 and  7 the  most  useful. 

It  is  expected  that  a stress  analyst  or  engineer  will  supplement 
the  solutions  presented  with  other  solutions  commonly  needed  in  his 
specific  work  and  that  a tabulation  of  toughness  data  for  the  materials 
will  be  available.  The  Damage  Tolerant  Design  Handbook  listed  in  the 
references  is  the  best  single  source  available  for  fracture  data  for 
metals . 

In  keeping  with  the  above  aims,  information  was  gathered  from  many 
sources,  including  the  book  by  Wilhelm,  the  handbook  by  Tada,  the  book 
on  three-dimensional  problems  by  Parmeter,  ASTM  STP  381,  the  seven 
volume  series  edited  by  Liebowitz,  and  various  technical  papers  by 
many  other  eminent  scientists  too  numerous  to  individually  acknowledge 
here.  It  is  hoped,  though,  that  due  credit  has  been  given  at  the 
appropriate  places  in  the  text. 
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Chapter  1.  INTRODUCTION 


1.1  History 

As  early  as  1920,  A.  A.  Griffith  [1]  published  his  concept  of  crack 
propagation.  The  concept  stated  that  an  existing  crack  will  propagate 
if,  in  so  doing,  the  total  energy  of  the  system  is  lowered.  The  stress 
analysis  used  to  calculate  the  stored  elastic  energy  was  based  on  a 1913 
solution  by  Inglis  [2]  for  an  elliptical  hole  in  an  infinite  plane  under 
uniaxial  tension.  Griffith,  in  his  work  on  glass,  postulated  the  exist- 
ence of  flaws  prior  to  loading,  which  were  large  compared  to  atomic  or 
molecular  distances.  In  1945  Orowan  [3]  noted  extensive  plastic  deforma- 
tion on  fracture  surfaces  of  materials  which  had  failed  in  a brittle 
fashion.  Then,  in  1948,  Irwin  [4]  pointed  out  that  the  Griffith  type 
energy  balance  must  be  between  the  strain  energy  stored  in  the  specimen 
and  the  surface  energy  plus  the  work  of  plastic  deformation.  The  same 
arguments  were  stated  by  Orowan  [5] , who  showed  that  the  Griffith  condi- 
tion, when  modified  to  account  for  plastic  work,  was  a necessary  and 
sufficient  condition  for  predicting  brittle  fracture. 

Then  in  1955,  Irwin  [6]  suggested  and  in  1957  [7]  he  showed  that 
the  energy  approach  is  equivalent  to  a stress  intensity  approach  according 
to  whych  fracture  occurs  when  a critical  stress  distribution,  character- 
istic of  the  material,  is  reached.  This  discovery,  shifting  fracture 
mechanics  from  the  energy  approach  to  the  stress  approach,  placed  frac- 
ture mechanics  on  a basis  more  familiar  to  the  practicing  engineer,  and 
was  largely  responsible  for  the  subsequent  development  and  application 
of  the  new  displine. 


The  need  for  increased  knowledge  of  brittle  fracture  mechanics  was 
clear.  The  spectacular  failures , by  breaking  in  half , of  the  liberty  ships 
in  the  1940's , the  repeated  brittle  fracture  failures  of  the  Polaris  missile 
at  stresses  below  the  design  level,  and  the  accelerating  demand  for  light- 
weight structures  provided  considerable  impetus  for  study  of  the  subject. 

By  1959  the  so-called  Griff ith-Irwin  concept  of  sharp  crack  "mechanics 
had  become  well  known  and  the  ASTM  Special  Committee  on  Fracture  Testing 
of  High  Strength  Materials  was  formed  to  launch  a broad  assault  on 
fracture.  By  then,  the  need  to  design,  and  test  specimens  with  severe 
artificial  flaws  was  recognized  and  there,  followed  much  experimental 
plane  strain  fracture  toughness  testing.  In  addition  to  the  experi- 
mental work,  considerable  analytical  work  was  performed,  seeking  to 
clarify  the  state  of  stress  around  flaws  of  various  shapes  under  various 
loading  conditions. 


In  summary,  the  history  of  fracture  mechanics  shows  development 
along  both  theoretical  and  experimental  lines.  From  a theoretical 
point  of  view,  the  subject  is  founded  upon  linear,  small  strain  elas- 
ticity insofar  as  the  stress  field  which  controls  the  crack  is  concerned. 
For  a given  material  crack  propagation  occurs  when  the  intensity  of  the 
stress  field  reaches  a certain  critical  value,  determined  experimentally. 
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1 .2  Philosophy  and  Purpose 
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The  basic,  underlying  philosophy  in  fracture  mechanics  lies  in  the 
assumption  that  the  elastic  stress  field  in  the  neighborhood  of  the 
crack  root  controls  the  crack  behavior.  The  influence  of  this  stress 
field  can  be  measured  in  terms  of  the  so-called  stress  intensity  factor 
(referred  to  as  K or  as  SIF),  or  alternatively,  from  a strain  energy 
approach,  in  terms  of  the  strain  energy  release  rate  g.  An  understanding 
of  of  these  quantities  and  their  correlation  is  basic  to  an  understanding 
of  this  approach.  In  addition,  it  is  assumed  that  of  all  flaws  in  a 
given  specimen  or  structural  component,  only  one  of  the  flaws  plays  a 
major  role  in  inducing  fracture.  In  other  words,  the  mechanism  requiring 
a coalesence  of  voids  into  a nucleus  out  of  which  a crack  is  created  is 
not  included  in  the  mathematical  modeling  which  follows.  This  is  done 
to  permit  comparison  with  controlled  experiments  on  precracked  materials 
and  also  because  fracture  of  a structural  component  is  complete  when 
only  one  crack  propagates  across  it. 


It  must  be  further  emphasized  that  the  subsequent  analysis  will  be 
accurate  only  for  brittle  fracture  resulting  from  sharp  cracks.  When 
applied  to  notches  with  root  radii  large  compared  to,  say,  the  grain  size, 
the  plastic  zone  preceding  the  crack  root  may  not  be  small  relative  to 
the  elastic  zone  which  is  influenced  by  the  notch  root.  In  such  cases, 
the  elastic  stress  field  surrounding  the  plastic  zone  may  not  govern; 
such  is  the  case  with  ductile  cracking.  However,  as  already  pointed  out, 
the  energy  of  plastic  work  is  not  negligible  in  technologically  important 
brittle  fracture  situations,  as  was  originally  assumed  by  Griffith. 


1 . 3 Energy  Rate  Analysis  of  Crack  Extension  and  Its  Measurement 

Before  proceeding  further,  it  seems  appropriate  to  present  an  ele- 
mentary analysis  together  with  a qualitative  explanation  of  a quantity 
called  the  "strain  energy  release  rate  "g,"  or  alternatively  the  "crack 
extension  force."  This  quantity  represents  a concept  which  may  be  new 
to  some  engineers,  and  although  the  interpretation  of  this  quantity  will 
become  more  sophisticated  in  later  chapters,  the  following  quasiheuristic 
elementary  discussion  will  serve  to  provide  the  basic  concepts  involved. 


Consider  a linearly  elastic  body  which  contains  a through  crack 
(here  defined  as  two  planar  mating  surfaces  which  are  stress  free  but 
are  initially  separated  by  an  infinitesimally  small  distance).  The 
body  is  in  a deformed  equilibrium  position  under  the  action  of  load  P, 
which  causes  an  extension  A under  the  load  point.  The  energy  rate  || 
(that  is,  energy  per  unit  of  new  crack  area  generated),  g,  available 
for  an  increment  of  crack  extension,  da,  is  provided  from  work  done  by 
the  force,  PdA,  and  the  release,  -dV,  in  the  total  strain  energy,  V, 
stored  in  the  body.  In  other  words,  part  of  the  work  PdA,  results  in 
an  increase  in  strain  energy  and  part  is  used  in  creating  new  fracture 
surface.  Expressed  mathematically 
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dA  dV 
da  da 


(1.1) 


8 = 


strain 
area  (■ 


energy  release  per  unit  increase  in  fracture  surface 

in. - lb 

. 2 ’ 
in. 


V = total  strain  energy  stored  in  the  body  (in. -lb)  as  indicated  in 
Figure  1.1.  (There  is  no  provision  in  this  explanation  to 
account  for  plasticity  but  this  is  a detail  to  be  dealt  with 
later. ) 


Figure  1.1.  Cracked  body  under  load. 


Measurement  of  g is  accomplished  by  the  definition  of  compliance, 
\ = /P,  an  inverse  spring  constant.  In  Equation  (1,1),  V = r PA 


-i>2' 


Thus , 


dV 


da  ~ P 5 


dP 


P2  b\ 
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then  by  substituting  into  Equation  (1.1) 


(1.2) 


>p 

Since  the  quantity  v-  cancels  from  Equations  (1.2),  g is  independent  of 

the  type  of  test  fixture,  whether  fixed  cross  head  or  fixed  load.  This 
fact  will  be  useful  later. 


In  a compliance  calibration,  a displacement  gage  is  mounted  on 
the  specimen.  A calibration  curve  of  X versus  a is  prepared  by  obtaining 
A - P curves  for  a specimen  with  various  crack  lengths  a (Figure  1.2). 


a = 0 a = 1 


Figure  1.2.  Compliance  calibration  curves. 

The  P - A curve  is  linear  and  its  recipocal  slope  is  the  X for  that 
particular  crack  length.  The  X - a curve,  however,  is  nonlinear.  For 
a specimen  with  a given  crack  length  or  area,  the  quantity  oX/oa  can  be 
determined  from  the  X - a curve.  Using  this  with  a known  load  P in 
Equation  (1.2),  the  strain  energy  release  rate  g becomes  known. 

Although  Griffith's  fracture  criterion  was  originally  stated  a 
slightly  different  way,  essentially,  crack  extension  occurs  when  g 
approaches  a critical  valve,  e.g.,  gc>  equal  to  the  total  solid  state 

surface  tension,  T.  That  is,  when 

g = gc  = T (1.3) 

then  crack  extension  occurs. 
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Instead  of  measuring  g as  explained  previously,  from  Griff ith's  work 
[1,8],  one  can  calculate  g for  the  plane  stress  condition  and  the  geometry 
shown.  Figure  1,3.  For  that  case. 


g = 


2 

2rtO  a 


where  E is  Young's  modulus  and  a is  the  applied  stress. 


(1.4) 


Figure  1.3.  Large  thin  plate  containing 
a crack  of  length  2a. 

An  incremental  crack  extension  of  da  on  each  end  of  the  crack  creates  a 
new  surface  area  of  4da ; hence  in  Equation  (1.3)  let  T = 4U,  where  U is 
the  surface  tension  per  unit  area.  Using  Equation  (1.3),  a gross  criti- 
cal (or  fracture)  stress  can  be  obtained. 


- =/— 

c v :ra 


(1.5) 
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It  is  noted  that  the  fracture  stress  is  inversely  proportional  to  the 
square  root  of  crack  length.  As  noted  before,  Irwin  and  Orowan  dis- 
covered that  for  most  materials  the  energy  balance  is  primarily  between 
the  elastic  energy  release  and  the  plastic  work  in  crack  propagation, 
and  have  suggested  the  addition  of  a plastic  work  factor  (which  is  large 
compared  to  U)  to  the  surface  tension  U. 

1.4  Experimental  Correlation 

It  has  already  been  mentioned  that  the  fracture  theory  is  based 
upon  linear  small  strain  elasticity.  This  is  feasible  if  Griffith's 
assumption  regarding  the  finiteness  of  the  crack  length  in  size  relative 
to  the  molecular  dimensions  holds.  It  further  permits  one  to  ignore 
atomic  inhomogeneities  in  formulation  of  the  theory.  Consequently, 
except  for  the  obvious  presence  of  a stress  singularity  at  the  tip  of 
a sharp  crack,  the  theory  will  be  in  many  ways  similar  to  those  for 
dealing  with  deep  notches. 
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This  means  that  experimental  models  which  must  be  tested  for  corre- 
lation of  theory  with  reality  will  be  much  the  same  as  models  one  might 
fabricate  for  choosing  theoretical  elastic  solutions  for  deep  grooves, 
(i.e.,  ordinary  continuum  theory  is  assumed  to  apply).  However,  there 
is  one  important  difference.  In  the  groove  problem,  no  stress  singularity 
exists,  and,  if  the  experiment  can  be  carefully  performed  for  one  material 
without  exceeding  the  elastic  range,  exact  correlation  with  theory  should 
result  and  it  may  be  assumed  to  hold  for  all  materials.  In  the  case  of 
a sharp  crack,  a plastic  zone  will  develop  due  to  the  stress  singularity 
in  the  experiment.  (One  can  only  insure  this  by  precracking  the  specimen, 
since  any  artificial  crack  may  not  have  the  appropriate  sharpness.) 

Since  the  plastic  behavior  is  different  for  different  materials,  tests 
must  be  run  for  each  material. 

Specimen  configuration  and  size  must  also  be  considered,  together 
with  the  loading  mode. 

1.5  Size  Effect 

Irwin  and  McClintock  [9]  have  identified  an  order  of  magnitude  size 
scale  appropriate  to  the  analysis  of  crack  behavior  phemomena  as  shown 
in  Figure  1.4  (all  dimensions  in  cm). 


Until  model  9 can  be  analyzed  by  means  of  wave  mechanics  and  the 
results  integrated  through  this  order  of  magnitude,  it  will  be  necessary 
to  consider  a series  of  different  models,  each  of  which  is  appropriate 
to  its  own  level  of  investigation  (i.e.,  size  range). 

Boundary  conditions  must  be  obtained  from  the  next  largest  and  next 
smallest  scales.  The  smaller  the  model  size,  the  broader  the  range  of 
conditions  which  can  be  predicted  or  correlated  from  given  experimental 
data.  From  whatever  level  one  starts,  test  data  must  be  tsed  in  place 
of  the  lack  of  knowledge  at  the  smaller  size  levels.  Note  that  the  local 
stress  field  surrounding  an  elastic  singularity  would  contain  on  the 
6 9 12 

order  of  (10°)  or  10  separate  dislocation  fields  where  one  field  is 
regarded  as  the  largest  size  which  would  be  sensitive  to  a single 
dislocation, 

1.6  Loading  Modes 

Stress  fields  near  crack  tips  can  be  divided  into  three  basic  types, 
each  of  which  is  associated  with  a Local  mode  of  deformation.  They  are 
known  as  Modes  I,  II,  and  III,  respectively  (Figure  1.5). 


I 


L* 


OPEN  MODE  TENSION 


EDGE  SLIDING 


EDGE  TEARING 


MODE  I 


MODE  II 


MODE  III 


Figure  1.5.  Modes  of  crack  tip  deformation. 


Modes  II  and  III  bear  a certain  analogy  to  edge  and  screw  disloca- 
tions respectively  in  the  sense  that  displacement  discontinuities  exist 
along  the  crack  surface  behind  the  crack  tips.  However,  the  components 
of  the  stress  vector  must  also  vanish  on  the  crack  surfaces  whereas  for 
dislocation  they  take  on  continuous  non-zero  values. 
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The  superposition  of  these  three  modes  is  sufficient  to  describe 
the  most  general  case  of  local  crack  tip  stress  and  deformation  fields. 
Furthermore,  it  is  clear  from  elasticity  theory  that  if  the  specimen 
width  B is  small  and  the  loads  are  distributed  uniformly  or  symmetrically 
over  B,  then  Modes  I and  II  would  essentially  conform  to  plane  stress 
fields  but  if  B is  large,  the  states  of  deformation  of  these  modes  might 
be  more  nearly  plane  strain.  Mode  III  is  clearly  neither  one. 

1.7  Summary 

From  the  foregoing  introduction,  the  following  points  may  be  made: 

a)  The  theoretical  foundation  for  fracture  mechanics  may  be 
soundly  based  upon  the  linear  small  strain  theory  of  elasticity,  regard- 
less of  the  phenomena  occurring  within  the  plastic  zone  at  the  crack 
tip  which  precipitates  fracture. 

b)  The  failure  criterion  may  be  based  upon  a limiting  intensity 
of  the  local  elastic  stress  field  in  the  neighborhood  of  the  crack  tip. 
This  limit  can  be  specified  in  terms  of  a single  material  parameter 
(not  yet  shown).  This  parameter  was  initially  recognized  historically 
to  be  the  so-called  strain  energy  release  rate  and  was  later  shown  to 
be  related  to  the  stress  intensity  factor. 

c)  The  correct  value  of  the  stress  intensity  factor  to  be  applied 
in  design  should  be  based  upon  appropriate  experimental  data  in  conjunc- 
tion with  the  appropriate  theoretical  stress  analysis. 
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Chapter  2.  THE  CRACK-TIP  STRESS  FIELDS 


2.1  Westergaard  Stress  Function 


The  most  direct  approach  to  the  determination  of  the  stress  and 
displacement  fields  associated  with  the  three  modes  of  loading  of  the 
crack  is  based  upon  the  use  of  the  Westergaard  Stress  Function  [10,11], 
Modes  I and  II  can  be  analyzed  as  plane  strain  problems  which  are  sym- 
metric and  skew-symmetric,  respectively,  with  respect  to  the  crack  plane. 
Mode  III  can  be  regarded  as  the  pure  shear  or  torsion  problem  (antiplane 
strain).  The  Westergaard  Stress  Function  will  be  defined  and  then  used 
in  defining  the  stress  intensity  factor.  Correlation  with  the  stress 
concentration  factor  and  the  energy  rate  analysis  will  follow. 


It  is  known  in 
$(x,y)  defined  by 

elasticity 

theory  [ 

12] 

1 that  any  stress  function 

_ ^2(t 

_ ;-2$ 

2. 
c $ 

Jxx  > 2 ’ °yy 

°y 

' ox2  ’ ** 

’ cxc-y 

(2.1) 

will  satisfy  the  equilibrium  equations 

of 

elasticity.  Furthermore, 

if 

the  function  satisfies  ~ i = 0 at  every  point  in  the  body  and  on  the 
boundary,  then  the  function  $ is  said  to  be  the  solution  to  the  associ- 
ated plane  strain  problem.  The  stress  components  are  shown  in  Figure  2.1. 


Figure  2.1,  Coordinates  measured  from  the  leading  edge 
of  the  crack  and  the  stress  components  in  the  neigh- 
borhood of  the  crack  tip. 


Now  by  differentiating  the  first  of  the  above  with  respect  to  x and  the 

? 

second  with  respect  to  y it  is  seen  that  v (Re  Z)  = 0.  Similarly, 

"7^  (1m  Z)  = 0.  The  above  equations  will  be  used  in  evaluating  the 


ji  derivatives  of  the  Z functions  in  the  following  paragraph. 


2.2  Mode  I - Stresses  and  Displacements 

For  Mode  I deformations,  in  order  to  free  the  crack  surface  of 

t , consider  the  following  special  form  of  Equations  (2.2). 
xy 

= Re  Zj.  + y Im  Z^  (2 

where  the  functions  Re  Z ^ and  Im  Z are  harmonic.  The  subscript  "I" 
refers  to  Mode  I.  To  get  the  stresses  use 


xx  o 2 

y 


o®  6 (Re  Z ) _ C(lm  Z ) 

^ — — + Im  ZT  + y — 

oy  oy  I oy 


— - Im  + Im  Zj.  4-  y Re 


= y Re  Z. 


— = Re  Zx  + y 


c(Re  Zx) 


= Re  Zj  - y Im  Z| 


Hence 


0 = Re  Z - y Im  Z'  . 

xx  I I 


Proceeding  similarly  one  obtains 


o = Re  ZT  - y Im  Z' 
xx  I I 


0 = Re  Z,  + y Im  Z' 

yy  I I 


t = -y  Re  Zl 
xy  I 

Thus  the  choice  of  the  form  of  3>  above  provides  that  t =0  every- 

I xy 

where  y = 0.  Furthermore,  using  Hooke's  Law  for  a state  of  plane 
stra in, 


T-  n.--  * 
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du  o v(o  + a ) 

x xx  yy  zz 


xx  ax 


ou  a v(o  + a ) 

x yy  xx  zz 


yy  °y 


! /OU  o u \ 0/1  , 

e = I _*  + ^ 

xy  2 l oy  ox  / E 


(2.7) 


e = 0 = 
zz 


a - v(a  + a ) 
zz  xx  yy 


where  u and  u are  the  displacements  in  the  x and  y directions,  respec- 
x y 

tively;  E is  Young's  modulus  and  v is  Poisson's  ratio. 

In  order  to  obtain  expressions  for  the  displacements,  substitute 
Equations  2.6  into  Equations  2.7  an  integrate  as  indicated.  For 
example 

ou  o - v(o  +CT) 
v _ yy xx  zz 

oy  E 


and  since  a = v(o  + ct  ) 
zz  xx  yy 


if  = E [(1  " y2)  'yy  ' V(1  + v)  °xxl  = [(1  - v)  CT: 


-VO 

yy  xx  J 


(1  - v)  (Re  ZI  + y Im  Zj)  - v (Re  Zj.  - y Im  Zpj 


1 + v 

E 


^(1  - 2 v)  Re  Z1  + y Im  Z|' 


Integrate  by  parts  and  use  Equation  (2.4) 


u = - g [(1  - 2 v)  Im  Zj.  + y (-Re  Zj.)  + / Re  ZJ  dy 


1 + v 
E 


(1  - 2 v)  Im  Zz  + y (-Re  Z^)  + Im  zJ 


uy  = ^ g V 2(1  - v)  Im  Zj.  - y Re  Zj. 


(2.8) 


S imilarly , 


u 

x 


1 + V 
E 


|(1  - 2v)  Re  Z. 


- y Im  Z 


(2.9) 


When  a straight  through-the-thickness  crack  lies  on  the  x axis  in  an 
infinite  plate  and  loads  are  applied  remote  from  the  crack  so  as  to 
produce  the  symmetric  Mode  I deformation  near  the  crack  tip,  one  may 
take 

<Dj.  = ReZj.+  ImZj 

Whence, 

o = Re  Z,  + y Im  Zl 

yy  I I 

t = - y Re  Z' 
xy  J I 


Thus,  t =0  along  y = 0.  One  must  also  require  that  Re  ZT  vanish  on 
’ xy  I 

the  crack  surface  in  order  to  make  a =0  there. 

yy 

In  addition,  note  that  (a)  stress  singularities  must  exist  at  the 
ends  of  the  crack,  and  (b)  u must  be  double  valued  along  the  crack 

surface.  When  a function  is  analytic  and  single-valued  in  a domain 
except  along  one  line  along  which  it  becomes  double-  or  multi-valued 
and  discontinuous,  that  line  is  called  a branch  cut.  If  one  takes 

z = 

l(z  + b)(z  - a)]1^ 

i 

then  for  a crack  extending  from  x = -b  to  x = +a  as  in  Figure  2.2  stress 
singularities  exist  at  x = -b  and  x = +a.  If  one  requires  that  g(z)  is 
well  behaved  (no  singularities,  etc.)  in  the  neighborhood  of  the  crack 
j tip,  then  the  crack  surface  can  be  throught  of  as  a branch  cut.  It  is 

also  required  that  Im  g(x)  = 0 for  -b  < x < a so  that  Re  Z = 0 in  that 

i range  (i.e.,  consider  as  a special  case  of  g(z)  = u(x,y)  + iv(x,y), 

| g = iz).  Then,  along  y = 0 for  -b  < x < a,  the  denominator  of  Z^  above 

is  imaginary;  also  g = ix  is  imaginary,  so  Z is  real  and  * 0.  Thus, 

■ 0.  So  one  must  have  Im  g(x)  = 0. 

4 yy  ° 

Moving  the  origin  to  the  crack  tip  at  x = a,  introduce  f = z - a 
(Figure  2.2).  Now 


Figure  2.2.  A crack  in 
a large  plate. 


7 = 8(C  + a» 

1/2 

1 Ur  + a + b)  M i/ 


which  can  be  written  as 


z = iiU 

ZI  ?l/2 

where  f(*)  is  also  well  behaved  so  that  when  t -»  0,  f(£)  remains  finite. 
Such  a well-behaved  analytic  function  may  be  expanded  into  a MacLaurin 
series  about  the  new  origin, 

00 

f (O  = f (0)  + E f(n)(0)  £n 
n=l 

write  the  results  as 

iUy2r-  + ^ + a2  (;  + . . . . 

ZI  = ^172 

where  f(0)  = (the^/Eit  is  arbitrary),  and  a^  is  f^(0).  Thus 


I:  ' ;^o  (2jx^)1  2 


(2.10) 


All  crack  tip  stress  functions  reduce  to  the  above  form.  The  functional 
forms  of  the  local  elastic  stress  and  displacement  fields  may  be  obtained 
by  substituting  the  above  into  the  stress  and  displacement  field 
equations. 

For  example,  for  it  is  seen  that  Re  Z^.  and  y Im  Z j.  are  required. 

Introduce  polar  coordinates  r and  / and  consider  the  origin  to  be  at  the 
right  crack  tip.  i.e., 


Substituting  into  the  expression  for  a 


Equations  (2.6) 


I , 9 1 . a . 39. 

~Yj2  (cos  2 " 2 Sln  Sln  ~2  ^ 


By  using  j sin  9 = sin  — cos  — , the  stress  0^  can  be  written  in  the 
desired  form.  Stresses  o and  follow  in  a similar  manner: 

yy  xy 

'S:  ‘ . 9 . 3 >. 

o = cos  — (1  - sin  r sin  — ) 

XX  rz 2 2 2 

y2nr 


*1  6 , . 9 . 39. 

a = ■—  cos  — (1  + sin  r sin  — ) 
yy  2 22 


I 9.9  39 

t = cos  r sin  7 cos  — 

Xy  Jin  2 2 2 


0 = v(0  + 0 ) T =T  =0 

zz  xx  yy  xz  yz 


where  is  the  stress  intensity  factor  for  Mode  I deformation.  The 

quantity  is  a constant  for  a given  loading  condition  and  a given 

specimen  geometry  and  crack  length;  it  will  be  discussed  in  more  detail 
later. 

Displacements  are  calculated  from  Equations  (2.8)  and  (2.9).  For 
example,  for  u^,  it  is  seen  that  Re  and  y Im  Z ^ are  required.  Obtain 

first  Re  ZT 


, . *1 
1 v^Ttl/2 


x/2^  J Jin 


- fl  „ 1/2  . 9 , . , 9. 

Zj  = / — Kj  r (cos  j + l sin  j) 


glUW.ik  .* 


Therefore, 


Re  ZT  =7-  Ki  r cos  - 


.1/2  _ 0 
2 


Now  1m  Zj  was  noted  during  the  calculation  of  the  stress  fields 


K 


Im  Zj  = - 


0 

1/2  SLn  2 


V'2n  r 

Use  y = r sin  9 to  obtain 


, „ KI  1/2.0., 

y Im  Z = - r sin  — sin  9 

x/2n 


Substituting  Re  and  y Im  Z^  into  Equation  (2.9)  yields 

.1/2  r 


1/  r 

= 2(1  + v)  I 


fix 


/i  o'!  8,1,9.- 

(1  - 2v)  cos  — + 2 sln  J Sln  ® 


The  displacement  u would  be  obtained  similarly.  In  order  to  write  u 

y x 

in  its  final  form  use 

sin  r?  = 2 sin  j cos  j 

and 


G = 


2(1  + v) 


Finally, 


K r 1/2 
KI  r -9 

u = cos  j 


x G 
1/2 


KI  r .9 
uy = sin  2 


y2ir  G 


, o^.20 

l-2v+sin  j 


2(1  - v)  - cos2  j- 


(2.12) 


u = 0 (for  plane  strain) 
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Equations  (2.12)  give  the  displacements  associated  with  Mode  1 
deformation.  Again,  it  is  pointed  out  that  different  Z^'s  change  only 

the  intensity  of  the  displacement  fields  so  long  as  first  mode  deforma- 
tion is  maintained. 

Note  also  that  the  obtained  values  of  u and  u involve  symmetric 

x y 

and  unsymmetric  forms  of  /2  and  so  they  do  indeed  describe  opening  mode 

deformations,  i.e.,  u Ct  sin  3/2  and  u a.  cos  0/2  and  so  are  unsymmetric 

y x 

and  symmetric,  respectively  (when  unsymmetric  forms  are  squared,  as 

? 

sin“  3/2,  the  unsymmetric  effects  vanish). 

2 . 3 Mode  II  - Stresses  and  Displacements 

Just  as  Equation  (2.5)  leads  to  the  proper  stress  and  displacement 
form  for  Mode  I deformation,  one  can,  by  a judicious  choice  of  4>  , obtain  the 

appropriate  forms  for  a Mode  II  deformation.  Once  $ has  been  chosen, 

the  procedure  followed  is  the  same  as  that  followed  in  the  development 
just  completed  for  Mode  I displacements.  Consequently,  instead  of 
developing  the  results  in  detail,  only  an  outline  of  procedure  with  the 
results  will  be  given.  Assume, 

®ri  = -y  Re  Zn  • (2.13) 

The  stresses  in  terms  of  become 


Mext,  taking 


uh r ‘0  (2rtf)1/2 


(2.16) 


to 

where  ” = re  , and  using  the  same  reasoning  as  before,  one  obtains 


Si  . • ■ 39 , 

°xx  = ' T/2  sin  2 (2  + COS  2 C°S 


y/2n  i 

Si 


yy 


y/2ri  r 
K. 

Xy 


30 

YJ2  cos  2 sln  ~ cos  ~ 


II  9 ,,  . 0 . 30. 

~J2  cos  J (1  - sin  j srn  — ) 


(2.17) 


v(c  + a ) , t = x =0  (plane  strain) 
xx  yy  ’ xz  yz  ^ 


and 

K -1/2 

II  . ....  2 0. 

u = sm  r (2  - 2y  + cos  r-) 

X x/2?  G 2 2 

„ 1/2 

^11  r 29 

u = cos  - (1  - 2v  + sin  -)  (2.18) 

> ^/2k  G 

u = 0 (plane  strain) 


where  K is  the  stress  intensity  factor  for  Mode  II  deformation. 

Equations  (2.18)  show  that  u is  an  unsymmetric  function  of  0/2  and  u is 

x y 

symmetric  function  of  9/2,  which  is  in  accordance  with  the  Mode  II  deformation. 

Due  to  the  forms  exhibited  by  and  4>  , Equations  (2.5)  and  (2.13), 

the  results  from  the  first  and  second  modes  may  be  directly  superimposed 
since  i>j.  + 5 is  also  an  admissable  form  of  Equation  (2.2).  It  follows 

that  the  stresses  and  displacements  resulting  from  individual  values  of 
and  ® may  also  be  directly  superimposed. 
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2.4  Mode  III  - Stresses  and  Displacements 

The  antiplane  strain  problem  or  two-dimensional  problem  of  pure 
shear  may  be  specified  as  follows: 

u = 0 
x 


u = 0 

y 


(2.19) 


uz  = uz(x,y) 


Applying  the  strain-displacement  equations  and  Hooke's  Law,  one  finds 


eu  ou  t 
2 _ z + __  x _ xz 

xz  cx  oz  G 


(2.20) 


CU  OU  T 

ey  oz  G 


so  that 


t = G - — 
xz  ox 


t = G -r— 
yz  oy 


(2.21) 


All  other  components  of  strain  and  stress  vanish,  and  the  single  result- 
ing equilibrium  equation  is 


= 0 

-X  oy 


(2.22) 


Substituting  Equations  (2.21)  into  Equation  (2.22),  there  results 


-2  - n 

v u =0 
z 


(2.23) 


The  function 


u = - Im  ZTTr 
z G III 


(2.24) 


where  , being  analytic  and  harmonic,  satisfies  the  governing  equation, 

Equation  (2.23).  Now,  following  the  same  reasoning  as  before,  the 


function  Z ^ for  a crack  along  the  x axis  must  take  the  following  form 
near  the  crack  tip: 


•o  (2*  n 


where  is  the  Mode  III  stress  intensity  factor.  Now,  substituting 

Equation  (2.25)  into  Equations  (2.21),  where  * = re  , and  proceeding 
as  before  one  obtains 


1/2  Sln  2 


1/2  COS  2 


a = c = a = t 

xx  yy  zz  xy 


(1  1/2  . 

- r sin  j 


ux  = u = 0 . (2.27) 

Having  now  obtained  the  expressions  for  the  stresses  and  displacements 
in  terms  of  the  stress  intensity  factors  , K^,  for  Modes  I,  II, 

and  III  loading,  it  is  pointed  out  that  one  is  able  to  do  this  only 
because  the  crack  surfaces  are  the  dominating  influence  on  the  stress 
distributions  in  the  neighborhood  of  the  cracks,  and  the  specimen  bound- 
aries and  loads  affect  only  the  intensity  and  not  the  form  of  the  local 
elastic  stress  and  displacement  fields. 

The  stress  Equations  (2. 11)  and  displacement  Equations  (2.12)  for  the 
Mode  I deformation  and  the  stress  Equations  (2.17)  and  displacement 
Equations  (2.18)  for  Mode  II  deformation  have  been  written  for  the  case 
of  plane  strains  (i.e.,  u = 0)  but  by  setting  a =-0  and  replacing 

2 Z Z 

Poisson's  ratio  v withv/(l+v)  (G  is  not  changed)  in  the  displacement 
equations,  field  equations  for  plane  stress  may  be  obtained. 
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2 . 5 Characteristics  of  the  Elastic  Crack-Tip  Stress  Field 

The  dominant  mathematical  feature  of  the  three  stress  fields  is  the 

1/9 

1/r  singularity.  All  three  sets  of  stresses  (Modes  I,  II,  and  III) 
were  obtained  by  neglecting  the  higher  order  terms  in  r in  the  MacLaurin 
expansion  of  Z;  thus,  the  stresses  can  be  regarded  as  good  approximations 
when  r is  small  compared  to  the  other  planar  dimensions  (x  - y plane)  of 
the  body,  such  as  crack  length,  and  mathematically  exact  as  r approaches 
1/2 

zero.  The  1/r  singularity  predicts  that  the  stresses  approach 
infinity  as  r approaches  zero.  Crack-tip  plasticity  prevents  this  in 
a way  which  is  discussed  in  Chapter  4.  Equations  (2.11),  (2.17),  and 
(2.26)  also  predict  that  the  stresses  approach  zero  as  r becomes  large. 
This  is  because  the  higher  order  terms  in  the  MacLaurin  expansion  were 
omitted.  Physically,  for  large  r,  the  stresses  would  approach  those 
values  in  the  body  that  would  exist  due  to  remote  loads  even  if  the 
crack  were  not  present. 

It  is  instructive  to  take  a more  detailed  look  at  the  nature  of  the 
crack-tip  stress  field  represented  by  Equations  (2.11)  for  the  opening 
mode.  The  various  stress  quantities  have  been  summarized  by  Williams  [13] 
as  shown  in  Figure  2.3.  It  is  seen  that  a state  of  hydrostatic  tension 
exists  in  the  material  along  the  crack  plane  (y  = 0)  ahead  of  the  crack 
tip:  i.e., 


a (x,0) 
xx 

a (x,0) 

yy 


= l 


(2.28) 


This  indicates  a decreased  tendency  for  yielding  at  the  base  of  the 
crack.  At  the  same  time  a large  amount  of  distortion  is  indicated  off 
to  the  sides  of  the  crack  (1  = ±70  degrees).  The  maximum  shear  stress 
occurs  at  ‘ = 90  degrees.  The  maximum  tensile  stress  occurs  not  ahead 
of  the  crack  as  one  might  expect  but  at  9 = 60  degrees. 


The  stress  components  near  the  crack  tip  in  polar  coordinates 
(Figure  2.4)  can  be  calculated  by  substitution  of  Equation  (2.11)  into 
stress  transformation  expressions  such  as 


o_=o  sin 

XX 


2 

A+o  cos  9 - 2t  sin  9 cos  9 

yy  xy 
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Figure  2.4.  Stress  components  in  polar 
coord inates . 

For  the  opening  mode,  Figure  2.5  shows  a sketch  of  the  variation  of  these 
stresses  with  the  angle  9.  The  location  of  maximum  cr  is  important  in 

determining  the  direction  of  crack  propagation,  and  Yoffee  [14]  assumed 
propagation  would  occur  along  a line  normal  to  maximum  a . For  the 

opening  mode,  maximum  ± occurs  at  9 = 0 degrees.  For  a mixture  of 

Modes  I and  II,  maximum  c would  occur  at  r 0 degrees. 


Figure  2.5.  Polar  plot  of  stresses. 

2 . 6 Meaning  of  the  Stress  Intensity  Factor 

The  primary  parameter  in  fracture  mechanics  is  the  stress  intensity 
factor;  it  is  important  to  understand  its  meaning.  For  a given  Wester- 
gaard  stress  function  Z note  that  from  Equation  (2.10), 


VL,  = lim  (2s  *)X  2 Z. 

^ I r 1 


(2.29) 


with  similar  equations  for  and  K . It  is  pointed  out  that  while 

will  be  different  for  different  Z^'s,  the  form  of  the  stresses. 

Equation  (2.11),  for  any  Mode  I problem  will  be  the  same.  The  only 
result  will  be  that  the  intensity  or  magnitude  of  the  stresses  in  the 
neighborhood  of  the  crack  tip  will  be  altered.  In  other  words,  KT  is 

l/2  1 

the  strength  of  the  1 r stress  singularity,  and  it  is  not  a function 

of  coordinates  r and  9. 

It  is  important  to  distinguish  here  between  the  stress  concentration 
factor  and  the  stress  intensity  factor.  The  stress  concentration  factor 
is  defined  as  the  maximum  stress  at  the  tip  of,  say,  a notch  divided 
by  the  nominal  stress.  Since  the  opening  mode  stresses  of  Equations  (2.11) 
mathematical  approach  infinity  at  the  crack  tip,  the  stress  concentration 
factor  for  a crack  would  be  infinite  while  the  stress  intensity  factor  is 
a finite  constant  (for  a given  specimen  geometry  and  load ) . This  constant  is 
multiplied  times  the  form  of  the  stress  field.  While  the  stress  concen- 
tration and  stress  intensity  factors  are  different  quantities,  mathema- 
tically, they  are  related  in  a useful  way  (Art.  2.8). 


l 


V •■>*** 


Finally,  from  dimensional 
(2.11),  (2.17),  and  (2.26),  it 
vary  linearly  with  the  applied 
also  depend  upon  the  geometric 

The  dimensions  on  K^.,  K^,  and 


cons iderations  of  the  stress  Equations 
is  clear  that  the  stress  intensity  factors 
loads  (for  linearly  elastic  bodies)  and 
configuration,  including  the  crack  size. 

1/2 

K are  (stress)  (length) 


An  aspect  of  the  stress  intensity  factor  which  has  caused  consider- 
able confusion  must  be  explained.  As  a holdover  from  the  formulative 
period  of  fracture  mechanics,  there  exists  another  definition  of  K - 
call  it  K*  - which  differs  from  the  present  definition  by  the  factor 

y/'n.  That  is,  K = x/ n K*.  When  the  K*  definition  is  used  then  the  singu- 
larity of  the  stress  fields,  Equations  (2.11),  (2.17),  and  (2.26)  are 

written  in  the  form  — rather  than  — . Unfortunately,  many  authors 
\'2r  v/^c7 

do  not  explain  which  definition  is  being  used,  and  one  must  examine  the 
form  of  the  author's  stress  field  for  clarification.  To  avoid  an  error 

of ,/x"when  using  stress  intensity  factor  solutions  from  the  literature 
one  must  be  careful  to  determine  which  definition  the  author  used. 


2 . 7 Dimensional  Considerations  for  the  Classical  Problems  of  Cracks 
in  Infinite  Plates 


In  Figure  2.6(a)  only  two  characteristic  dimensions  are  present, 
and  a.  Furthermore,  this  configuration  is  symmetric  with  respect  to 
the  crack  plane.  Due  to  this  symmetry  with  respect  to  the  x-z  plane, 
one  notes  that  only  Mode  I deformation  is  present.  This  leads  to 

1/2 

Equation  (2.111  for  stress.  Since,  in  Equation  (2.11),  a K/r  , 
one  sees  from  dimensional  considerations  that 


KI  = C1  ° a 


1/2 


(2.30) 


Also,  due  to  symmetrv,  K. 


II 


"i 


II 


0.  For  this  case  the  complete 


Westergaard  stress  function  is  available  so  that  the  exact  stresses 
ahead  of  the  crack  can  be  computed.  The  exact  stresses  can  then  be 
compared  with  the  factor  mechanics  stress  Equations  (2.11)  to  determine 
the  constant  c^.  For  the  case  of  part  (a),  the  Westergaard  stress 


tunction  is, 


Z = 


oz 


, 2 2,1/2 

(z  - a ) 


(2.31) 
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Figure  2.6.  Crack  in  an  infinite  plate  - 
various  modes. 


From  Equation  (2.6), 


ReZj  + ylmZj.  . 


Along  y = 0 


i = Re  Z = 


yy i„_r 


I 

(J  + 


+ a) 

2 2 
- 4 r - a 


where  ( is  a complex  coordinate  centered  at  the  right  crack  tip. 
£ = re1  and  for  y = 0,  i = 0 and  £ = r.  Therefore, 


a(r  + a) 


o(r  + a) 


=0  T,  ,2  2]1  2 [2  l1 

(r  + a)  - a J (r  + 2ar) 


is  the  exact  stress.  From  Equations  (2.11),  the  fracture  mechanics 
stress  is 


=0  (2nr) 


By  requiring  that  the  approximate  stress,  Equation  (2.33)  equals  the 
exact  stress,  Equation  (2.32),  as  r -*  0 


r -»  0 (2itr ) 


there  results 


K = a (;ta) 


o(r  + a) 

, „ ,1/2  1/2 

(r  + 2a)  r 


(2.34) 


Comparing  Equation  (2.34)  with  Equation  (2.30),  it  is  seen  that 

C = /r 

2.8  Correlation  of  Stress  Intensity  Factor  with  Stress  Concentration 
Factor 

Many  solutions  are  already  available  in  the  literature  which  provide 
stress  concentration  type  solutions  for  various  kinds  of  notches.  It  is 
then  a distinct  advantage  to  be  able  to  relate  the  stress  concentration 
solutions  for  a notch  to  the  stress  intensity  factor  for  a crack. 

Consider  a symmetrically  loaded  notch  for  which  the  tips  are  embedded 
in  a Mode  I stress  field.  For  this  case,  it  is  known  from  Equation  (2.34) 
1/2 

that  Kj.  = c(rta)  . From  Equation  (2.11)  the  stress  a directly  ahead 

1/2 

of  a crack  (9  = 0)  is  c = Kj./(2itr)  . Suppose  one  used  this  form  to 

approximate  the  stress  near  a notch  with  a root  radius  p and  root  stress 
(this  is  accurate  only  when  p ->  0).  Then  in  order  to  have  the  proper 

units  on  J one  must  pick 

yy  r 


K a c' 


(2.35) 


Now  by  expressing  J in  terms  of  a from  the  notch  solution  and  letting 

■ 0,  in  the  limit  one  obtains  an  expression  for  in  terms  of  0 and 

a.  By  equating  this  new  expression  to  a known  case  Equation  (2.34),  c' 
can  be  evaluated.  Consider  the  notch  solution  for  an  elliptical  hole 
under  uniform  tension. 


o = c 1 + 2 (-) 

m c 


(2. 36' 


where  p,  the  semiminor  axis,  becomes  the  root  radius  as  p becomes  small. 
Using  Equation  (2.35)  as  p becomes  small 
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K = 1 im 

p -*  0 


'c'  cr 


1 + 2 [- 

P 


1/21 


1/2  I 


= 2 c ' oa 
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Equating  this  to  Equation  (2.34),  it  is  seen  that 

c .-4 


or  that  Kj  for  a crack  under  Mode  I loading  can  be  computed  from  a 
corresponding  notch  solution  by 


*1 


= linu 


4 a. 


(2.37) 


A converse  relationship  also  exists  for  determining  stress  concentration 
factors  for  notches  from  stress  intensity  factor  solutions  for  cracks 
[15]. 

A similar  line  of  reasoning  may  be  used  to  identify  the  problems  of 
Figure  2.6(b)  and  (c)  with  Modes  II  and  III,  respectively,  so  that: 

Figure  2.6(b),  Mode  II: 

Kjj.  = T(*a)1/2 

Kx  = Kii  = ° . 


(2.38) 


Figure  2.6(c),  Mode  III: 
Km  = T(*a)1/2 


K 


•'ll  " ° 


(2.39) 
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2.9  Equivalence  of  Energy  Rate  and  Stress  Intensity  Factor  Approaches 

It  was  shown  in  Chapter  1 that  g is  independent  of  the  load  type 
or  grip  condition.  Thus,  for  convenience,  in  the  discussion  which 
follows,  the  fixed  grip  situation  may  be  employed  with  no  loss  in 
generality  of  results. 

If  an  elastic  body  is  loaded  and  the  grips  (load  point  displacements) 
are  then  fixed,  the  strain  energy  change  dV/da  will  be  the  only  contri- 
bution to  g.  In  this  case,  the  work  required  to  close  a small  segment 
of  the  crack  5a,  from  the  opened  position  (Figure  2.7(a))  to  the  closed 
position  (Figure  2.7(b))  is  identical  to  the  strain  energy  released 
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Figure  2.7.  Closing  of  a crack  by 
crack-surface  tractions. 

upon  its  opening  (for  linear,  elastic  behavior).  The  work  can  be  com- 
puted as  the  crack-surface  tractions  required  to  close  the  crack  times 
their  closing  displacements  times  one-half.  The  tractions  will  be  those 
along  a for  case  (b)  of  Figure  2.7.  The  displacements  will  be  those 
along  5a  in  case  (a)  of  Figure  2.7. 


The  -train  energy  release  per  unit  area  of  crack  extension  is: 

6a 

2_ 
i 

0 


dV 

g = ' dl 


= lim  I 

fixed  5a  -»  0 J 

grips 


O U XU  T u 

yy  y . *y  * . yz  z 
2 2 2 


dx  . 
(2.40) 


The  negative  sign  is  properly  accounted  for  because  the  crack  opening 
displacements  which  are  used  here  are  negative  to  the  crack  closing 
displacements.  The  factor  of  2 preceding  the  integral  is  necessary 
because  displacements  occur  on  both  sides  of  the  crack  plane.  The 
stresses  for  Modes  I,  II,  and  III,  Equations  (2.11),  (2.17),  and  (2.26), 
are  evaluated  at  0 and  r = x.  The  displacements  for  the  three 

modes  are  evaluated  at  = x and  r = 5a  - x.  Working  out  the  details, 
the  strain  energy  release  rate  becomes, 


S 


(2.41) 
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Equation  (2.41)  indicates  a separate  energy  rate  contribution  from  each 
of  the  three  stress  field  modes.  Then  write, 


8 8I  + 8II  + 8III 


where 


- 1 ' V v2 
8II  E *11 

„ = 1 + u J- 

6 III  E III 


(2.42) 


(2.43) 


Equations  (2.43)  has  been  developed  for  a state  of  plane  strain.  They 

2 

can  be  converted  to  the  plane  stress  case  by  discarding  the  1 - U 
term  in  the  first  two  equations. 

Thus,  for  isotropic  materials.  Equations  (2.43)  permit  the  conver- 
sion of  a measured  g into  a K without  knowledge  of  the  stress  function 
for  the  problem  at  hand.  It  was  assumed  in  obtaining  Equations  (2.43) 
that  the  crack  extends  straight  ahead  along  a plane.  For  combined-mode 
loading,  cracks  do  not  extend  in  a plane  fashion  and  therefore  results 
obtained  using  the  approach  taken  to  obtain  Equations  (2.43)  would  be 
divergent  from  reality.  Equations  (2.43)  do  serve,  though,  to  show  the 
general  equivalency  of  the  energy  rate  and  stress  intensity  factor 
approach. 

It  was  discussed  briefly  in  Chapter  1 that  fracture  is  assumed  to 
occur  when  the  strain  energy  release  rate  g approaches  a critical  value 
which  depends  upon  the  material;  i.e.,  when  g -*  It  can  now  be  seen 

that  this  is  equivalent  to  assuming  that  fracture  occurs  when  the  stress 
intensity  factor  approaches  a critical  value.  Then  the  fracture  criterion 
is  that  as 

K - Kc  , (2.44) 

crack  propagation  occurs. 
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Chapter  3.  STRESS  INTENSITY  FACTOR  SOLUTIONS 


The  relationship  of  the  stress  intensity  factor  to  the  Westergaard 
Stress  Function  was  described  in  Chapter  2.  The  discussion  will  begin 
here  bv  considering  some  examples  of  solutions  by  the  Westergaard  method. 
Its  limitations  will  be  noted,  and  a more  general  approach,  based  upon 
the  method  of  Muskhelishvili  will  be  developed  for  plane  problems. 
Finally,  a brief  discussion  will  be  included,  noting  the  results  of  other 
technologically  important  solutions. 

3 . 1 Examples  of  Crack-Tip  Solutions  Using  the  Westergaard  Stress 
Function 


Consider  a crack  of  length  2a  in  a uniform  biaxial  stress  field 
(Figure  3.1).  This  problem  was  treated  from  a dimension  viewpoint  in 
Chapter  2.  The  Westergaard  Stress  Function  for  this  case  is  given  by 
[10,111 
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Figure  3.1.  Crack  in  a uniform,  biaxial 
stress  field. 


(3.1) 
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Use  now  z = £ + a to  write 
_ 2 o (C  + a) 


1 " R + 2„l/V/2 

and  use  Equation  (2.29)  to  yield 

K = lim  (2nC)1/2  ° ^ \ n Yil 

IChO  a + 2a)w  C 

There  results 


Kj  = 0(ia)l,‘ 

which  was  found  previously  in  Equation  (2.34). 

Consider  next  an  infinite  periodic  array  of  cracks  along  a line 
in  a field  of  uniform  tension  (Figure  3.2).  The  half  period  is  b and 
the  half  crack  length  is  a.  The  stress  function  that  solves  this  pro- 
blem has  been  provided  by  Westergaard  [10,11]  as: 


(sing)  - (sing) 


231/2 


Figure  3.2.  Periodic  array  of  cracks  along  a line. 


Again  using  Equation  (2.29)  with  z = £ + a one  finds 


, .1/2  /2b  :ta\ 

“ c(na) 


This  formula,  sometimes  referred  to  as  the  tangent  solution,  has 
been  widely  used  in  fracture  toughness  testing.  It  is  a good  example 
of  how  solutions  for  infinite  bodies  can  sometimes  be  usefully  applied 
to  finite  shapes.  Consider  the  following. 

The  dotted  lines  on  the  plate  in  Figure  3.2  represent  axes  of 
symmetry  and  so  they  are  free  of  shear  stress.  Now  if  a « b,  then  one 
may  remove  the  horizontal  stress  field  at  infinity  leaving  only  small 
self -equi librating  stresses  a along  these  axes.  This  will  not  alter 

the  value  of  since  depends  only  upon  a . Consequently,  one  may 

regard  it  as  permissable  to  cut  the  sheet  along  the  axes  of  symmetry 
and  to  consider  the  above  as  an  approximate  solution  for  finite 

width  strips  with  central  cracks  provided  a < b/2.  Results  of 
Isida  [16]  and  Kobayashi  [17,18]  indicate  that,  even  for  a > b/2,  the 
error  in  this  practice  is  less  than  7 percent. 

Furthermore,  if  one  slices  the  body  along  the  y-axis  (or  correspond- 
ing axes)  the  above  may  be  regarded  as  an  approximate  solution  for 

double-edge-notched  strips  which  is  within  2 percent  according  to 
Bowie  [ 19]  if  a > b/2. 

Both  the  center-cracked  and  double  edge-cracked  specimens  have  been 
used  in  fracture  toughness  testing.  Finally,  note  that  for  b ».a  the 
results  found  here,  Equation  (3.4),  reduce  to  the  case  of  a single 
crack  in  an  infinite  plate.  Equation  (3.2). 

3.2  Crack-Tip  Solutions  Using  the  Method  of  Muskheli shvi li 


Westergaard  stress  functions  have  been  used  to  solve  many  problems, 
but  there  are  limitations  to  the  scope  of  the  method.  The  method  is 
restricted  to  infinite  two-dimensional  bodies  with  cracks  along  one 
or  more  straight  lines.  The  complex  function  approach  developed  by 
Muskhelishvili  [20]  and  others  [11]  has  some  advantages  over  the 
Westergaard  method  by  treating  a broader  class  of  plane  problems. 

As  noted  in  Chapter  2,  an  Airy  stress  function  $ must  satisfy 

: = 0 together  with  the  boundary  conditions  of  the  problem.  It  is 
known  from  elasticity  theory  (see  Timoshenko  and  Goodier  [12])  that  $ 
can  be  written  as 

* = Re[  zi(z)  + \(z)]  (3.5) 


where  T(z)  and  X(z)  are  suitably  chosen  complex  functions  known  as 

Goursat  functions  and z is  the  conjugate  of  the  complex  variable  z.  It 

is  further  known  that  the  stress  invariant  a + a is  given  by 

xx  yy 

a +a  = 4 Re[ J ' (z)]  (3.6) 

xx  yy 

where  '(z)  represents  dT(z)/dz.  Introduce  now  the  concept  of  the 
complex  stress  intensity  factor 

K = Kj  - i K1Z  . (3.7) 


Now  adding  the  first  two  of  Equations  (2.11)  to  the  first  two  of 
Equations  (2.17)  one  obtains 


2K 


a + a 

xx  yy 


(2Jtr) 

Now  consider  the  function 


T72  cosf 


2K 


(2*r) 


% Sin  2 


(3.8) 


f~T  FT  ( a . J a \ 
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It  is  seen  that 

T 


Re 


' = /^(Ki  COS  | - Kn  Sinf 


(3.10) 

which  is  the  same  as  the  right-hand  side  of  Equation  (3.8).  Therefore, 


a + a = Re 
xx  yy 


(Kt) 
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Now  let 


= z - z , 


(3.11) 


(3.12) 


where  is  the  location  of  the  crack  tip  (Figure  3.3).  Now  substituting 
Equation  (3.11)  into  Equation  (3.10)  one  obtains 
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0 +0 
xx  yy 


Re 


yir  (z  - z^ 
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Figure  3.3.  Complex  coordinates 
at  the  crack  tip. 

But  irora  Equation  (3.6) 

J +o  = 4 Re.  '(z) 

xx  yv 

Now  it  is  known  that  when  the  real  parts  of  two  analytic  complex  func- 
tions are  the  same,  then  the  functions  themselves  can  differ  by  only 
a real  constant  times  i.  It  is  also  known  from  elasticity  theory  that 
in  a singly  connected  body,  the  form  of  the  complex  stress  function 
for  a unique  stress  state  can  only  be  specified  to  within  a real 
constant  times  i.  For  unique  displacements,  even  this  constant  must 
vanish.  By  restricting  consideration  to  the  neighborhood  of  the  crack 
tip  (i.e.,  a singly  connected  region)  one  may  write,  without  loss  of 
the  above  prescribed  extent  of  uniqueness. 


= 4 1 ' ( z ) 


->  Zj  - =LV 


Therefore , 


2 v 2jt  lim  I 


(z)  • (z  - z^) 


(3.13) 


The  technique  oi  conformal  mapping,  to  map  holes  into  cracks,  has 
been  used  by  Sih,  Paris,  and  Erodogan  [21],  Muskhelishvili  1 20 ] , and 
others  to  determine  the  function  ,(z)  for  an  extensive  variety  of  crack 
problems . 

Consider  the  mapping  functions: 


z = -(  ,) 


(3.14) 


v..v  * 1 


where  >,  is  the  complex  variable  in  the  mapped  plane  and  the  usual 
restrictions  as  to  analyticity  and  single  valuedness  apply.  Then 


*•(«)  = = 

d’i  dz  W (l) 

and  Equation  (3.13)  becomes 

K = 2 lim  |:o(’|)  ” 


1/2  .'.Lil 

V ( i) 


Now,  corresponding  to  a crack  tip  at  in  the  z-plane,  there  will  be 

a crack  tip  at  ■,  in  the  , plane.  In  order  to  transform  a straight 

crack  of  length  2a  in  the  x-direction,  located  centrally  at  the  origin 
in  the  z-plane,  into  a circular  hole  of  unit  radius  in  the  rj  plane,  con- 
sider the  mapping  function: 


I 


or  in  terms  of  mapped  variables  u and  v (Figure  3.4), 


a ...  1 

iv  = - u + tv  + — 

2 u + iv 
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Figure  3.4.  Mapping  a crack  into  a unit  circle. 


Separating  Equation  (3.17)  into  real  and  imaginary  parts,  one  obtains 


a . u 

x = 2 U 2 + 2 
U V 


a v 

y - 2 v " 2 2 

u -f  v 


The  boundary  y = 0 is  transformed  as  follows.  Let  y = 0 in  the  second 
of  Equations  (3.18). 
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Solutions  are: 

v = 0,  a line  of  v =0,  and 
2 2 

u + v -1=0, 

a circle  of  unit  radius,  as  shown  in  Figure  3.4. 

To  locate  the  right  crack  tip  in  the  mapped  plane  let  z -»  z^  = a 
in  Equation  (3.16).  Then, 


solving 


i.e.,  u = 1,  and  v = 0.  Similar  reasoning  shows  that  the  left  crack 
tip  z = -a  maps  into  the  point  q^  = 1 or  u = -1,  and  v = 0. 


Now  using  the  mapping  function.  Equations  (3.16)  in  (3.15), 


Letting,  ^ = 1 and  taking  the  limit  as  q -»  q , one  finally  finds 

K = 2 /j  ; '(1)  . (3.19) 

Thus,  it  is  noted  that  stress  intensity  factors  may  be  found  simply  from 
a knowledge  of  ' ( ,)  in  the  vicinity  of  the  crack  tip.  While  only 
• ( |)  is  used  in  Equation  (3.19),  it  should  be  pointed  out  that  in  the 
course  of  solving  a new  problem  both  T(q)  and  X(q)  would  be  needed  in 
satisfying  the  boundary  conditions. 


As  an  application  of  Equation  3.19  consider  the  problem  shown  in 
Figure  3.5.  For  this  case  Muskhelishvili  [20]  has  given  the  stress 
function  as 


Figure  3.5.  Skewed  crack  in  a uniaxial  tensile  field. 
Taking  the  derivative 


Evaluating  at  t = 1 

*'(1)  = ~ (1  - ei2r) 
or 

ca  t 

(1)  - ~ (1  “ cos2p  - lSmr) 

■4 

Using  Equation  (3.19)  together  with  the  identity  1 - cos  2p  = 2sin2. 
and  separating  into  real  and  imaginary  parts,  one  finds. 


Kj.  j = a s/ r.a.  sin  cos  . 
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Another  important  example  of  a problem  solved  by  the  method  of 
Muskhelishvili  is  that  of  a straight  crack  with  a concentrated  load 
as  shown  in  Figure  3.6.  For  this  case,  Sih,  Paris,  and  Erdogan  [21] 
found  the  stress  function  to  be 


♦(,) 


< 

1 + K 


log  •,  - log  ( ,0 


(3.21) 


v 


Figure  3.6.  Straight  crack 
with  a concentrated  force. 


where 


F = P - iQ 

Ig  corresponds  to  z = b 

K = an  elastic  constant  3 - 4r  for  plane  strain  and  (3  - t>)/(l  + t>) 
for  plane  stress  (v  = Poisson's  ratio). 


Applying 


K 


I 


K 


II 


Equation  (3.19),  one  obtains, 


(3.22) 


Equations  (3.22)  are  important  because  they  become  the  Green  influence 
functions  for  solving  any  single  straight  crack  problem  in  an  infinite 
plane  from  a knowledge  of  the  stresses  on  the  prospective  crack  surface 
with  the  crack  absent;  i.e.,  a (x,o)  and  t (x,o).  The  procedure  is  to 
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solve  the  infinite  plate  problem  without  a crack  for  the  above  stresses 
at  the  prospective  crack  surface  location.  Then  superimpose  tractions 
equal  and  opposite  to  those  computed  on  the  crack  surface  location.  The 
result  can  be  written  directly  from  Equations  (3.22)  as: 


KI  = 


II 


(ita) 


1/2 


+a 

/ 


a (x,o) 

yy 


+a 


(-a) 


1/2 


/ 


(x,o) 


-a 


a + x 
a - x 


a + x 


1/2 


dx 


1/2 


dx 


(3.23) 


Equations  (3.23)  give  the  stress  intensity  factors  at  the  right  crack 
tip.  The  multiplier  of  1/2  disappeared  from  the  results  of  Equations 
(3.22)  because  in  this  case  loads  are  applied  to  both  crack  surfaces. 
The  shearing  contribution  due  to  Q or  (x,o)  disappeared  from  the 


first  equation  because  the  force  was  applied  to  both  surfaces 
for  the  contribution  of  P or  a (x,o)  in  the  second  equation. 

3.3  Influence  of  Local  Boundaries  Other  than  Crack  Surfaces 


similarly 


Thus  far,  only  cases  where  the  crack  location  was  r^vs^re  from  the 
external  boundaries  of  the  body  have  been  considered  so  that  infinite 
body  solutions  were  acceptable.  One  may  separate  problems  where  boun- 
daries other  than  the  crack  surface  boundaries  are  also  local  into  two 
categories: 


a)  Edge  cracks  in  semi-infinite  bodies. 

b)  Cracks  in  infinite  strips. 

Edge  Cracks  in  Semi-Infinite  Bodies  - It  was  implied  in  the  pre- 
vious article  that  good  anproximations  could  be  obtained  to  problems 
involving  cracks  in  semi-infinite  bodies  (b/a  » 1)  by  slicing  infinite 
body  problems  along  axes  of  symmetry.  The  accuracy  of  such  solutions  has 
been  evaluated  by  comparison  with  analytical  solutions  by  Bowie  [ 19] , 
Wigglesworth  [22],  and  Irwin  [23,24].  Although  dimensional  analysis  of 
.he  problem  in  Figure  3.7  would  lead  one  again  to  Equation  (2.9), 

1 '2 

j witn 

required  to  obtain  the  effects  of  the  free  surface  of  the  semi-infinite 
body.  In  general,  the  method  of  Muskhelishvili  is  employed,  in  which 
the  mapping  function  takes  the  form  of  a series  and  the  use  of  Green's 
influence  function  lead  to  dual  integral  equations.  Since  these  methods 
have  already  been  developed  and  illustrated  by  means  of  simple  problems, 
the  results  only  will  be  written  here  for  the  problem  of  Figure  3.7: 


Kt  — C,  a ( a)  *■  with  C.  unknown,  formidable  analytical  methods  are 
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: = 1.12  J(r.a) 


(3.24) 


Figure  3.7.  An  edge  crack  under 
Mode  I loading. 


Recalling  Equation  (3.2)  K = o(ira)  which  was  developed  for  an  inte 

nal  crack  in  an  infinite  body,  one  may  conclude  that  the  free  surface 
correction  factor  for  edge  notches  perpendicular  to  a field  of  uniform 
tension  is  1.12.  On  the  other  hand,  note  the  results  of  Sih  [25]  for 
an  edge  crack  (Figure  3.8)  under  a mode  III  deformation;  his  results 


are  identical  with  those  for  an  internal  crack  in  an  infinite  body  and 
therefore  no  free  surface  correction  is  required  for  this  case. 


Figure  3.8.  An  edge  crack  under 
Mode  III  loading. 


Since  the  free  surface  would  normally  be  a loaded  surface  in  the 
infinite  plate  with  Mode  II  loading,  it  is  concluded  that  there  is  no 
Mode  II  case  directly  analogous  to  the  Mode  II  case  in  the  infinite 
plate  corresponding  to  the  cases  discussed  above. 

Cracks  in  Infinite  Strips  - This  class  of  problems,  including  plate 
strips  with  transverse  internal,  edge,  and  dual  colinear  edge  cracks 
under  normal  tension  has  been,  together  with  the  same  configurations 
under  plane  bending,  of  great  practical  interest  in  connection  with 
fracture  toughness  testing.  As  is  often  the  case,  however,  closed  form 
solutions  for  such  problems  are  not  available  and  many  of  the  approxi- 
mate solutions  in  the  literature  may  be  of  questionable  accuracy.  Con- 
sequently, it  seems  important  not  only  to  site  available  results  but  to 
give  estimates  of  their  accuracy  as  well. 


The  limitations  on  Equation  (3.4) 


= o(jta) 


/ 2b  Tta\ 

\Ta  t3n  2b/ 


(for  a row  of  colinear  cracks  in  an  infinite  plate)  when  it  is  applied 
to  cracks  in  strips  was  discussed  in  the  first  of  the  chapter.  These 
limitations  were  compared  with  solutions  resulting  from  direct  attack 
upon  strip  type  problems. 

One  of  the  most  useful  approaches  to  this  class  of  problems 
appears  to  be  the  work  of  Isida  [16].  Isida  has  extensively  developed 
mapping  functions  for  strip  problems  for  the  determination  of  stress 
concentrations  at  the  tips  of  round-ended  cracks  of  end  radius  p.  His 
results  are  presented  in  the  form  [16], 


(3.26) 


where  a = crack  half-length,  and  b = strip  half-width,  f is  deter 

mined  by  applying  power  series  mapping  to  the  stress  function  variable 
z and  : is  itself  expressed  as  a power  series. 


Equation  (3.26)  ma\  be  substituted  directly  into  Equation  (2.37), 


K = Um 
I p -*  o 


o p 

max 


* ..'ta.'V.S#  *•  a 


I 
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The  resulting  stress  intensity  factors  can  be  computed  to  any  degree  of 
accuracy  b Isida's  methods,  provided  the  power  series  employed  in  the 
analysis  converge,  which  they  do  for  relatively  large  variations  in 
a/b.  Within  this  minor  limitation,  Isida's  results  are  expected  to 
yield  accuracies  of  ±2  percent.  Idisa  l 16]  has  computed  results  in  the 
form  of  Equation  (3.26)  for  a number  of  problems  of  special  interest 
in  fracture  testing.  Consider,  for  example,  the  problem  of  Figure  3.9. 


C± 

a 

b 

a 

b 

* * 
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Figure  3.9.  Center  cracked 
infinite  strip. 


SubstitULing  Equation  (3.26)  into  Equation  (2.37). 

r 

lim 

I — 

? 


1/2  1/2 
- 2ca  f 


O -1/: 


, ,1/2  , 

K = u(;ra)  f 


1/2 


(?) 


■(/) 


becomes  the  free  surface  correction  factor  for  the 


finite  width  strip,  the  approximate  form  for  which  is  given  by 
Equation  (3.4).  The  two  solutions  are  compared  in  Table  3.1. 


(3.27) 


Frequent ly  in  design  situations,  approximate  solutions  must  be 
employed,  and  sometimes  solutions  for  a certain  class  of  problems  must 
be  adapted  to  cases  for  which  they  do  not  exactly  apply.  The  com- 
parison in  Table  3.1  gives  an  idea  of  the  kind  of  approximations  one 
can  make  for  a given  amount  of  error.  For  example,  for  a/b  as  large  as 
0.410  the  approximate  tangent  formula  gives  results  to  within  5 percent. 
Differences  for  larger  a/b  values  result  from  the  fact  that  lines  of 
symmetry  in  the  infinite  plate  are  subject  to  larger  local  stress  dis- 
tributions. 
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TABLE  '}.  U FREE  EDGE  CORRECTION 
FACTORS  FOR  CENTER  CRACKED 
INFINI  TE^TRI  PS 


0.074 
0.207 
0.275 
0.337 
C -.10 
0 . 46 o 
0.533 
0.592 


f(?) 

1 

.00 

1 

.03 

1 

.05 

1 

.09 

1 

.13 

1 

.18 

i 

.25 

1. 

.33 

/2b  na  V' 
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Baechner  [29]  has  developed  integral  equation  procedures  and  solved 
a number  of  crack  problems  including  the  solution  to  a strip  with  a 
single  edge  notcn  subjected  to  pure  bending  (Figure  2.10). 


Figure  3 . iu.  Strip  with  single 
edge  crack  subjected  to  pure 
bending. 


His  results  are  given  b. 


-gm 


6M  / a \ 

where  c = and  f^—  Its  given  in  Table  3.2.  One  sees  a large 

(h  - a) 

variation  in  f(a/h)  for  small  a/h.  For  larger  values  of  (a/h),  f(a/h) 
tends  to  become  constant,  approaching  the  value  it  would  have  for  a 
deep  notch,  | 30 J . 


TABLE  3.2.  STRESS  INTENSITY  FACTORS  FOR  A SINGLE-EDGE 
CRACKED  STRIP  LOADED  IN  BANDING 


a/h 

0.05 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60  (and  larger) 

ii) 

0.36 

0.49 

0.60 

0 . 66 

0.69 

0.72 

0.73 

Bowie  developed  polynomial  mapping  functions  for  use  with  the  com- 
plex stress  function  technique  to  solve  plane  problems,  such  as  cracks 
emanating  from  circular  holes  [31]  and  the  double-edge-cracked  strip 
in  tension  [19] . The  latter  example  (Figure  3.11)  may  be  compared  with 
the  tangent  formula,  Equation  (3.4).  This  comparison  may  be  achieved 
by  introducing  a correction  factor,  f (a/b),  into  Equation  (3.3).  Thus 
one  has 


^ = c(na) 


l/2 


/ 2b_  «a\  l/i 

y Tta  *"  n 2b/ 


(3.29) 


II 


= K 


III 


0 


Figure  3.11.  Double-edge-cracked 
strip  loaded  in  tension. 
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The  results  are  shown  in  Table  3.3.  From  this  table,  one  sees  that,  for 
low  a/b  values,  the  correction  factor  of  1.12  for  a crack  emanating  from 
a free  surface,  as  noted  in  Equation  (3.24),  is  present.  As  a/b 
increases,  its  effect  disappears  and  the  correction  factor  becomes  1.00. 
The  values  in  the  last  column  represent  the  correction  factors  which 
would  have  to  be  applied  to  the  solution  for  a single  internal  crack 
in  an  infinite  plate  to  account  for  the  two  deep  cracks  in  a finite  width 
strip . 


TABLE  3.3.  STRESS  INTENSITY  SOLUTION 
FOR  A DOUBLF-EDGE-C RAC KED  STRIP 
LOADED  IN  TENSION 


a/b 

‘(f) 

L/b  = 1.00 

‘(f) 

(same 

L/b  = l.OCuas 

te-#/2‘(f) 

L/b  — > oo 

0.1 

1.13 

1.12 

1.12 

0.2 

1.13 

1.11 

1.12 

0.3 

1.19 

1.09 

1.13 

0.4 

1.16 

1.06 

1.14 

0.5 

1.14 

1 no 

1.15 

0.6 

1.10 

1.01 

1.22 

0.7 

1.02 

1.00 

1.34 

0.8 

1.01 

1.00 

1.57 

0.9 

1.00 

1.00 

2.09 

shut  s t;  at  wnen  using  Equation  (3.2-+), 
for  a/b  < 0.5  and  Equation  (3.4),  K,  = 


2b  nal 

fe  C3n  2bJ 


The  foregoing 

1/2 

K = 1.12  c(  : a) 1 , 

1/2 

0(  'a)  , for  a/b  > 0.5  for  the  double-edge-cracked  problem,  the  error 

will  be  less  than  3 percent  for  L/b  > 3.  Such  information  can  fre- 
quently be  of  help  in  estimating  the  stress  intensity  factor  for  design 

problems . 


Numerical  methods  involving  the  boundary  collocation  of  the  stress 
function  have  been  developed  by  Kobayashi  [18]  and  Gross.  Srawley  [32] 
and  b\  Wilson  [33]  and  bv  Newman  [34],  to  name  a few. 
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Kobayashi  118]  created  the  infinite  strip  with  a central  crack 
using  Muskheli shvil i ' s complex  stress  functions,  and  collocating  equally 
spaced  points  in  the  sides  and  ends  of  the  strip. 


Gross  [ 321  solved  the  single-edge-cracked  strip  using  the 
Williams  ( 13]  eigen  function  representation  of  the  Airy  stress  function. 
His  results  are  given  in  the  usual  form 


1 /2 

Kj  — tr  ( a)  f 


(?) 


(3.301 


KII  ~ KIII 


where  the  single  edge  cracked  strip  is  shown  in  Figure  3.12.  Gross'  [ )2] 
numerical  results  are  compared  with  Bowie's  [19]  solution  for  the  double- 
and  single-edge-cracked  strip  in  Table  3.4.  This  table  shows  how  the 
stress  intensity  factor  is  increased  due  to  bending  in  the  single-edge- 
cracked  strip.  When  the  double-edge-cracked  strip  is  cracked  half  way 
across  (a/b  = 0.5),  the  correction  factor  is  approximately  1.15,  whereas 
when  the  single-edge-cracked  strip  is  cracked  halfway  across  (a/b  = 1.0), 
the  correction  factor  is  2.82  — over  twice  as  great. 


Figure  3.12.  The  single-edge 
cracked  strip  loaded  in 
tension . 
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TABLE  3.4.  CORRECTION  FACTORS  FOR  THE  SINCLE-EDGE- 
C RACKED  STRIP  COMPARED  WITH  THE  DOUBLE-EDGE- 
CRACKED  STRIP 


a/b 

a 

2b 

Si ng  Le-Edge 
Crac  .ed 
f(  a/b) 

C ross  [32] 

Double -Edge -Cracked 

fe  Lan  if | f<*  b> 

Bowie  [19] 

Single- Edge -Or acked 
f ( a/b ) 

Bowie  [19] 

0.10 

VTT 

1.1* 

1.12 

1.15 

0.20 

o.lo 

1.19 

1.12 

1.20 

0 . 30 

0 L3 

i .29 

1. 13 

1.29 

0 . 40 

0.20 

1.3' 

1.14 

1.37 

0.50 

0.25 

1.50 

1 . 15 

1.51 

0.60 

0.30 

1.66 

1.22 

1.68 

0.70 

0.35 

1.87 

1.34 

1.89 

0.80 

0.40 

2.12 

1.57 

2.14 

0.90 

0.45 
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2.09 

2.46 

1.00 

0.50 

2.82 

- 

2.86 

3 . -»  Three-Dimensional  Considerations 

The  discussion  thus  far  has  only  dealt  with  the  cwo-dimensior.al 
case.  In  practical  anolications  of  results  though  it  should  he  remem- 
bered that  all  bodies  are  really  three  dimensional.  It  is  not  unusual 
to  treat  certain  problems  as  two  dimensional,  analytically,  and  bv 
correlating  with  experiment,  to  assess  the  three-dimensional  effects  ir. 
the  real  problem.  The  implication  here  is  that  the  three-dimensional 
effects  are  alwavs  present,  and  are  yet  often  justifiably  negligible. 

0:  er  roole.is  unden ia'  1 tall  for  a tnree-c xmensional  appro at  Sj- 
is  tiie  case  of,  sav,  an  eliipitcal  flaw  embedded  in  a bod;.. 

Consider  the  stresses  in  the  neighborhood  of  the  curved  crac. 
rvrder  shown  in  Figure  3.1  Kassir  anc  Sih  [35jhave  si, own  that  ii  the 
stresses  are  expressed  in  terms  of  a set  of  polar  coordinates  r and 
defined  in  a plane  normal  to  the  edge  >f  the  crac'-  that  the  angular 
distri nutionand  universe  square-root  stress  singular i tv  of  the  plane 
roblem  (Equations  2.11,  2.17,  and  2.26)  are  preserved.  This  means 
tit  the  local  elastic  stress  field  can  be  described  in  terms  of  the 
tutee  stress  intensity  factors,  and  that  the  concepts  already  discussed 
for  two  dimensional  problems  appl\  here  also.  In  the  presence  of  all 
three  modes  the  following  fracture  criterion  in  three  dimensions  was 
postulated  by  Kassir  and  Sih  [35] 
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Figure  3.13.  Stresses  and  coordinates  in 
a plane  normal  to  a crack  edge. 


This  equation  shows  that  the  crack  propagation  will  occur  when  some  com- 
bination of  Kj,  Kjj  and  K reaches  a critical  value.  In  many  three 

dimensional  problems;  however,  only  Mode  1 may  be  present. 

Using  a method  employing  Fourier  Transforms,  Snaddon  [ 36J  solved 
the  problem  of  a circular  "penn\  shaped"  internal  crack  loaded  with 
an  internal  pressure  on  the  crack  boundaries  (Figure  3.14).  The  result- 
ing stresses  and  displacements  were  given  as  sets  of  integral  equations. 
Upon  superimposing  a dilitational  field  to  negate  the  pressure  distri- 
bution in  the  crac';, 

' zz  ’ 'J  rr 


The  solution  for  a circular  crack  under  opening  mode  tension  was 
obtained.  For  this  case,  it  was  found  that: 


The  penny- shaped  crack. 


The  two-dimensional  plate  solution  is  modified  by  the  factor  — . 

Sadouskv  and  Sternburg  [37]  and  Green  and  Sneddon  [38]  have  discussed 
the  problem  of  stresses  in  an  infinite  body  near  an  elliptical  cavity. 
When  the  ca\itv  is  degenerated  into  a crack,  certain  difficulties  arise 
in  the  computation  of  the  stresses  at  the  edge  of  the  crack.  Irwin  [39] 
subsequently  computed  the  stress  intensity  factor  at  any  location  on  the 
crack  border,  described  by  the  angle  g (Figure  3.15)  by  comparing  the 
displacements  obtained  by  Green  and  Sneddon  with  Equations  (2.12).  His 
resulLs  showed  that' 


-L3) 

£(.0 


Ku  ■ hn  ' 0 


lln'5+  (b) 


2C 
cos  g 


(3.32) 


where 


2 2 

[ 1 - (k  ) sin“  ] d 


s the  complete  elliptic  integral  of  the  second  kina  and 


l . b 

k = 1 - — 


Figure  3.15.  Elliptical  crack. 


Now  for  very  large b (i.e. , b *3 
Jt/2 

E(lO  = J cos  d =1 
o 


and  i f 


= :'/ 2 , Kt  = a(rta) 


1/2 


which  is  identical  with  Equation  (2.34)  for  a through  crack  of  length 
2a  in  an  infinite  plate.  In  effect,  the  above  may  be  thought  of  as 
extending  the  thickness  of  an  infinite  plate  to  infinity.  On  the  other 
hand,  if  one  sets  b = a,  then 


S -foo:.)1'2 


wii^  ,s  the  sane  as  Equation  (3.31)  for  the  circular  cracK. 

The  companion  problem  of  the  above,  the  case  of  a flat  elliptical 
crac  . in  an  infinite  medium  loaded  with  an  arbitrary  shear  load 
(Figure  3. 16),  was  solved  by  Kassir  and  Sih  [35].  The  solution  for  a 
stress-free  crack  under  uniform  loads  in  an  infinite  medium  can  ce 
obtained  as  the  superposition  of  the  stresses  for  uniform  loads  at 
infinity  on  a uncracked  medium  and  the  stresses  due  to  the  negative 
of  the  uniform  loads  applied  directly  to  the  crack  surface.  Conse- 
q entlv , for  purposes  of  calculating  the  form  of  the  singular  stresses 
it  suffices  to  apply  the  loads  in  question  directly  to  the  crack  surface, 
as  shown  in  Figure  3.16. 
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Figure  3.16.  Elliptical  crick  under 
uni  lorn  shear. 


The  stresses  x and  x (Figure  3.13)  contained  only  the  tearing 
nt  tz 

mode  influence,  but  all  other  components  exhibited  a combination 

of  and  . In  general,  the  k's  were  a complicated  function  of  '0,<l>,u;, 

and  a/b . 

Semi-Elliptic  Surface  Flaw  — The  surface  flaw  is  important  for 
two  reasons:  (a)  it  has  caused  a number  of  well -documented  service 

failures,  and  (b)  it  has  been  used  extensively  in  fracture  toughness 
testing.  In  service  toe  surface  flaw  occurs  by  fatigue  growth  of  such 
things  as  a tool  scratch,  or  an  arc  strike.  This  flaw,  illustrated 
in  Figure  J.l7,  was  first  analyzed  b Irwin  [39].  By  modifying  the 
elliptic  flaw  solution,  Equation  (3.32),  Irwin  suggested  that  at  the 
deepest  point  of  the  flaw 


1.1  a 


w 


21 1/2 


(3.33) 


where  is  the  yield  strength  of  the  material.  At  the  deepest  point 

. 2» 

of  the  flaw,  = 90,  it  is  noted  that  che  bracketed  factor  in 
Equation  (3.32)  is  1.  The  negative  factor  in  Equation  (3.33)  is  a 
correction  for  crack  tip  plasticity  (See  Chapter  4),  and  the  factor 
of  1.1  is  an  approximate  correction  for  the  bisecting  free  surface 
which  transforms  the  elliptic  flaw  into  a semi-elliptic  flaw.  The 
factor  of  1.1  is  approximately  equal  to  the  edge  crack  correction  fac- 
tor of  1.12  previously  presented.  Since  Equation  (3.33)  contains  no 
correction  for  the  free  back  surface  of  the  plate  the  range  of  validit 
of  Equation  (3.33)  is  restricted  to  a/T  < 0.5. 


figure  3.17.  Semi-eliiptic  surface  flaw. 

Recentlv,  numerical  techniques  have  been  applied  ir.  order  to 
obtain  more  refined  solution  to  the  surface  flaw  problem,  tome  of 
these  results  are  given  ir.  Chapter  . Because  of  its  importance,  the 
surface  crack  was  the  subject  of  a special  symposium  of  the  19"?  Annual 
Meeting  of  ASME  [ uO] . 

Most  two-dimensional  problems  exhibit  three-dimensional  effects. 
Consider  the  problem  of  a sheet  of  finite  thickness  with  a through  crack 
(Figure  3.18).  For  an  infinitelv  thick  sheet,  a plane  strain  analvsis 
would  apply.  For  a vanl shingly  thin  sheet,  a plane  stress  analvsis 
would  anoly.  But  with  a finite  thickness,  a mixed  situation  of  plane 
stress  near  the  surface  of  the  plate  and  plain  strain  in  the  interior 
occurs  in  the  crack  tip  stress  field.  Consequently,  the  stress  inten- 
sity factors  computed  for  plane  problems  represent  only  values  averaged 
through  the  plate  thickness.  Recent  results  by  Leven  [41]  and  Cruse 
and  Van  Buren  [42]  reveal  a thin  boundary  laver  on  the  plate  surfaces 


. .■***>< 


-V*  - V*—  s 


} PLANE  STRESS 


I PLANE  STRAIN 
| PLANE  STRESS 


Figure  .n.  Thiriness  • an  at  ion  •>£  a through 

e r dc 


• Figure-  3.1.)  Ip  which  o drops  iron  the  '.air  strain  value  of 

2 z 

•(c  + » to  zero.  Accompanying  fhis  is  a sharp  drop  in  the  in-plane 

stresses,  and  . Consequently . tue  stress  intensity  factor  is 

X" 

grater  in  f’-e  interior  of  the  plate  than  on  the  surface. 

Altho  i crack  tip  plasticity  further  complicates  the  situation, 
the  above  indicates  in  part  why  a crack  begins  to  grow  in  the  interior 
of  a plate  rather  than  at  the  surface,  forming  a "tongue"  or  "thumb 
nail"  shaped  crack  front.  Even  though  this  effect  is  frequently 
observed  in  practice,  and  in  fact,  has  been  directly  utilized  in  the 
so  called  "pop-in"  tests,  ignoring  this  effect  will  not  lead  to  sig- 
nificant inaccuracies  of  computed  stress  intensity  factors  for  use  in 
developing  fracture  criteria. 

3 Estimation  of  the  Stress  Intensity  Factor 


In  design  a plications  a solution  ma\  not  be  available  for  the 
r den  at  and.  With  some  practice,  however,  it  becomes  possible  to 
i st i ate  . c stress  .c  isit  factor  for  cases  of  practical  interest. 

Consider  the  f llowirtg  example  from  Haris  and  Sih  [11].  The  crack  has 
an  arbitral-,  share,  si  own  in  Figure  3.19,  and  is  subjected  to  a uniform 
• •»!••*’  rci!  -.c  the  ac  "l.nt  Separate  regions  of  the  crack 

M . - , N,,  I".  . For  C is  CU  ■ » R__  = 0.  '.sing  tie 

r \ .ilc  tor  'e  ienn\ -shaded  crac  , 


and  tie  tunnel  crac 
1/2 


K 


r ai 


I 


V • 


one  can  guess  the  Mode  I stress  intensity  factor  for  the  different 
labeled  portions  of  the  crack  border.  Guess,  e.g., 


1/2 

k « 0.75  aCita^  ' 
\n 

k2  ~ 0.85  a(na^)  L 

k3  *>  1.05  0(rta3)1/2 
1 /2 

k,  = 0.75  a(ita  ) ' ^ 
4 4 


Ihese  estimates,  it  is  believed,  are  accurate  to  within  10  percent  and 
probably  even  5 percent. 


preceding  page  b 
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Chapter  4.  EFFECTS  OF  CRACK-TIP  PLASTICITV 

Up  to  this  point,  the  entire  effort  has  been  devoted  to  the  estab- 
lishment and  use  of  a linear  elastic  theory  for  the  local  stress  field 
around  the  tip  of  the  crack.  The  discussion  will  start  here  bv  con- 
sidering the  general  features  of  small  scale  plastic  yielding  at  the 
crack  tap.  Two  analytical  models  will  be  reviewed  and  their  results 
compared.  Some  remarks  will  be  made  concerning  the  limitations  of  the 
crack-tip  stress  theory. 

4 . 1 Introduction 

The  local  elastic  stresses  predicted  by  Equations  12.11),  (2.17), 
and  (2.26)  approach  infinity  at  tne  tip  of  the  crack.  For  a real 
material  it  is  known  that  the  stress  does  not  become  infinite  This 
means  that  yielding  must  occur  which  eliminates  the  crack-tip  stress 
singularity.  In  the  study  of  brittle  fracture,  this  yielding  is  what 
is  normally  referred  to  as  small  scale,  meaning  that  the  crack-tip 
plastic  zone  is  small  compared  to  the  elastic  stress  field  within  which 
it  is  embedded.  The  size  and  shape  of  this  plastic  zone  has  been 
studied  analytically  and  experimentally  by  several  researchers.  Due 
to  the  complex  analysis  required  for  Modes  I and  II,  an  exact  solution 
has  been  obtained  for  the  Mode  III  case  only. 

For  the  Mode  I case,  an  estimate  can  be  made,  following  McClintock 
and  Irwin  [9].  Figure  4.1  shows  the  locus,  R,  of  points  at  which  the 
Mises  yield  criterion  would  be  exceeded  if  the  elastic  stress  distribu- 
tion were  unaffected  by  yielding,  assuming  Poisson's  ratio  to  be  1/3. 

The  difference  in  the  plane-strain  and  plane-stress  loci  is  not  due  to 
the  in-plane  stress  components,  whicn  are  identical,  but  to  the  presence 
or  absence  of  transverse  constraint,  It  is  remembered  that  for  plane 

stress  ■ = 0,  whereas  for  plane  strain,  o = 0(0  + a ).  The 

z z r z z xx  yy 

presence  of  o causes  a state  of  triaxial  stress  which  inhibits  plastic 
r zz  r 

flow  for  the  plane-strain  case.  The  above  results  are  estimates  only 

since  yielding  arrects  tne  stress  distribution  outside  tne  plastic  zone 

and  thus  modifies  the  shape  of  the  plastic  zone  itself.  Material  work 

hardening  also  affects  the  shape  of  the  plastic  zone  and  the  extent  of 

yield ing . 


From  Figure  4.1,  it  is  noticed  that  the  radius  R of  the  plastic 

2 

zone  is  proportional  to  (K^/tT  : ) The  plastic  zone  size  is  normally 

stated  by  giving  its  radius  R^  ahead  of  the  crack  (i.e.,  at  4 = 0). 

For  this  dimension  the  commonly  assumed  values  are 


(plane  stress) 


Figure  '+.1  The  shape  of  the  Mode  I plastic  zone  determined 
from  miss  yield  criterion,  with  n = 1/3. 


R 


y 


6it 


(plane  strain) 


(4.1) 


where  is  the  yield  stress  of  the  material.  Since  the  yielding 

y s 

causes  a redistribution  of  stresses  ahead  of  the  crack,  the  actual 
diameter  is  estimated  at  twice  the  above,  i.e.,  2 Ry. 


Liu  [43]  proposed  a model  which  shows  the  plane-stress,  plane- 
strain  transition  of  a plastic  zone  in  a finite-thickness  plate.  The 
model,  indicating  the  relative  sizes  of  the  plastic  zones,  is  shown 

in  Figure  -*.2 


Based  on  the  work  of  Hahn  and  Rosenfield  [44],  Wilhelm  [45]  sum- 
marized the  effect  of  work  hardening  as  follows: 

a)  Large  work  hardening:  Dumbell  or  cirular  zones  are  formed 

!’rifjre  4 2' 

b)  Small  work  hardening:  Slip  is  concentrated  ahead  of  the 

crack  and  yielding  is  confined  to  a thin  wedge  ahead  of 
the  crack. 


Before  proceeding  to  a more  detailed  study  of  plastic  yielding, 
it  is  appropriate  to  note  how  a simple  plasticity  correction  can  be 
made  in  fracture  analysis.  It  will  be  shown  in  Article  4.2,  for 
Mode  III  yielding,  that  the  elastic  stress  field  in  the  presence  of  a 
plastic  zone  is  the  same  as  the  elastic  case  where  the  crack  tip  has 
been  moved  forward  by  half  the  width  of  the  plastic  zone.  Irwin  [ 46 ! 
proposed  using  an  effective  crack  length  determined  by  replacing  a 
with  a + Ry.  Then  expressions  for  K of  the  form 


0 


CRACK 

FRONT 


Figure  4.2.  Mode!  or  through-the-thickness  plastic 
zone  as  first  visualized  by  Liu 


K = o f (/I) 

dq  - jr.  .v 


K = ' f (v'a  + Ry) 

This  was  the  basis  of  the  plasticity  correction,  which  appeared  in  the 
surface  flaw  solution  Enuation  (3.33).  As  mentioned  previously,  the  com- 
monly used  values  of  Ry  are  those  given  by  Equat ions  (4. 1 ) . This  approach 
was  once  used  widely  in  fracture  toughness  testing,  but  has  now  been 
largely  replaced  by  the  use  of  compliance  techniques,  which  determine 
the  effective  crack  length.  In  fracture  analysis,  the  approximate  tech- 
nique of  replacing  "a"  with  a + Ry  is  still  used  although  for  many 
materials  and  situations  Ry  is  sufficiently  small  to  neglect  Indeed, 
for  high  accuracy  it  is  necessary  to  have  Ry  small  enough  to  neglect 
entirely.  For  low  accuracy  calculations  of  larger-scale  yield  problems, 
the  practice  of  replacing  "a"  with  a + Ry  is  regarded  as  a reasonable, 
approximate  method;  the  accuracy  depends  upon  the  extent  of  vielding, 
though  results  of  such  calculations  must  be  interpreted  with  judgement. 
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4.2 


Mode  III  — Elasto-Plastic  Solution 


Because  of  the  relative  simplicity  of  the  Mode  III  case  Rice  [47] 
was  able  to  find  a solution  to  the  plasticity  problem  for  a perfectly 
plastic  material . Although  th“  Mode  III  case  cannot  be  directly  adapted 
to  the  more  important  Mode  1 case,  the  results  have,  by  analogy,  been 
used  to  help  understand  certain  features  of  Mode  I plasticity.  Because 
of  its  importance,  the  solution  is  summarized  here 

Beginning  with  the  local  eaistic  stress  field  of  Eauations  (2.26), 
for  the  Mode  III  problem.  Figure  4 3(a), 
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(4.2) 


and  write 


T 

yz 


hn 

(x  - iy,1/2 


(4.3) 


s ince 


cos  y + i sin  y 


1 

- iy 


Now  seek  solutions  to  the  elasto-plastic  problem  by  specifying  that 
the  stress  aoproaches  the  value  of  Eauation  (4.3''  asymptotica  1 l-r  ir 
x - iy  — « . 


Define  a dimensionless  stress  by, 


t T 

yz  _ . _xz 
yp  ‘yp 


(4.4, 


where  is  the  vield  stress  in  shear 

yp 

can  write 


x - iy 


2 


Now  from  Equation  (4.3)  one 


(4.5) 
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Figure  4.3.  Small-scale  yielding  solution  for  elastic  - 
perfectly  plastic  Mode  III  deformation  (a)  physical 
plane  (b)  map  of  elastic  region  of  physical  plane  onto 
interior  of  unit  semicircle  in  stress  plane. 

One  .an  now  visualize  x - iy  as  a function  of  f.  The  asymptotic 
boundary  condition  can  then  be  stated  as  either 
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as 


x - IV 
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or 


x - iy 
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2x  T 
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as  i i ; 


(4.6) 


yp 


Ttv-  map  of  the  physical  plane  into  the  £ plane  is  shown  in  Figure  4.3. 
The  elastic-plastic  boundary  maps  into  a unit  semicircle.  Points  at 
infinity  map  into  the  origin  and  the  crack  surfaces  map  into  the 
imaginary  : axis. 


In  the  neighborhood  of  the  crack  tip  (i.e  , within  the  plastic 
zone),  the  stress  is  at  the  yield  level. 


f. 


xz 


yz 


yp 


or 


= 1 as  x - iy 1 -*  0 


(4.7) 


61 


1 _ .ASK*  7#**  ***&>  . 


.V*  v* 


Consider  the  function 

km  i,  . i } 
x - iy  = 1 1 + 2J 

P-l  - 5 

yp 


(4.8) 


t is  ,een  tnat  this  function  provides  the  singularity  at  infinity 
demanded  b the  boundarv  condition  Equation  (4  6)  and  leads  to 


near  the  crack  tip.  Th.  elastic-plastic  boundary  ray 

? i 

be  locate,  in  the  physical  plane  form  Equation  (4.3.'.  Since  - e 

on  the  elastic-plastic  boundary  in  the  stress  plane 


k; 

r.  - iy  --  — 1 1 + cos  2 - i sin  i -I 

2r.  t 

yp 


(4.9) 


which  is  the  path  in  the  physical  plane  corresponding  to  BED  in  the 
stress  plane  in  Figure  4.3.  Separating  the  real  and  imaginary  parts, 
one  has 


(1  + cos  2 • ) 
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yp 


(4.10) 


hll 
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2-  t 


yp 


These  are  the  parametric  equations  of  a circle.  Thus,  the  elastic- 
plastic  boundary  is  a circle  as  already  shown  in  Figure  4.3(a). 

The  diameter,  R , of  the  plastic  zone  is  given  by  x at  • = 0 in 
Equation  (4.10).  Thus, 


R = x 


w 


= 0 1 \ T 


yp 


(4.11) 


is  the  'lode  III  plastic  zone  diameter. 
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Solving  for  the  stresses  from  Equation  (4.8), 
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Theretore , 
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(4.12) 


are  the 
Equation 
a stress 

by  half 


stresses  in  the  elastic  region. 
(4.3),  it  can  be  seen  that  the 
field  identical  to  the  elastic 

the  plastic  zone  size  (i.e.,  by 


Comparing  Equation  (4.12)  with 
effect  of  yielding  is  to  create 
stress  field,  but  shifted  ahead 
o 2 

K Itt/'23X  t ),  as  if  the  crack 

7. 1 1 yp 


tip  were  at  the  center  of  the  plastic  zone  This  fact,  first  noted  by 
Irwin  [46]  and  Koskinen  [48],  lends  support  to  the  practice  in  Mode  I 
fracture  of  increasing  the  physical  crack  length  by  one-half  the  plas- 
tic zone  size.  Furthermore,  from  Equation  (4.11),  the  plastic  zone  size 
R depends  only  upon  KTTT  and  t . Thus  it  is  established  that  for 

O J r HI  yp 

sr.-.a  . 1 - sea  le  yielding,  the  plastic  zone  ; ize  uepenos  upon  tne 

elastic  stress  intensity  factor. 


4 . 3 Strip  Yield  Model 

In  lieu  of  an  exact  analysis  of  the  elasto-plastic  problem  for 
Mode  I loading,  a simple  model  was  proposed  by  Dugdale  [49].  This 
model,  called  a strip-yield  model,  provides  a known  plastic  zone  shape 
by  assuming  that  yielding  is  confined  to  a thin  layer  of  material 
directly  ahead  of  the  crack  tip.  Then  the  model  of  the  elastic-plastic 
problem  becomes  the  one  pictured  in  Figure  4.4.  While  a model  such  as 
this  connot  be  expected  to  yield  detailed  features  of  the  plastic 
deformation  near  the  crack  tip,  one  can  expect  reasonable  predictions  of 
such  things  as  plastic  zone  size,  functional  dependence  of  plasticity 
effects  upon  external  loads  and  geometry , and  unloading-reloading  behavior. 
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la/  EL  ASTO— PLASTIC  PROBLEM 


(bl  THE  STRIP- Y IEL  D MODEL 


STRESS 


W)  STRESS- STRAIN  CURVE  FOR 
RIGID-PLASTIC  MATERIAL 

Figure  4.4.  Strin-'-ield  modeling  of  the  crack-tip  plastic  zone. 

Consider,  as  an  example,  the  finite  crack  in  an  infinite  plate 

•mown  in  Figure  4.5.  In  that  figure,  the  portion  of  the  rigid  plastic 

material  which  has  yielded  has  been  removed  and  replaced  by  the  viald- 

stress  ' on  the  boundary  of  the  extended  portion  of  the  crack  w. 
vs 

Since  the  rigid  plastic  material  was  assumed  to  offer  no  resistance 
to  extension  or  contraction  in  the  x-direction,  the  solutions  to  the 
elastic  regions  of  the  model  in  Figure  4.5  may  be  obtained  by  caking 
the  usual  elastic  solutions  and  replacing  a by  a + w,  and  then  adding 
the  distribution  as  an  additional  external  loading. 

yp 
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PROTOTYPE  MODEL 

Figure  4.5.  Strip  yield  model  for  a finite  crack 
in  an  infinite  plate. 

The  problem  thus  formulated  can  be  solved  not  just  for  this  specific 
example  but  in  general  terms  which  will  cover  all  possible  Mode  I 
loadings  for  cracks  in  infinite  plates.  Let 

Kg(a)  = Elastic  stress  intensity  factor  for  the  crack  problem 

Ze(z,a)  = Westergaard  stress  function  for  the  crack  problem 

K (w,a)  = Elastic  stress  intensity  factor  for  the  a load 

P ys 


Zp (z ,w, a) = Westergaard  stress  function  for  the  s load 
where  the  subscript  I for  Mode  I has  been  dropped. 


Since  one  is  concerned  here  only  with  the  superposition  of  two  sets 
of  elastic  solution,  the  fact  that  the  physical  cracks  for  the  "p" 
solution  will  have  one  side  of  the  crack  running  into  the  other  will  be 
of  no  consequence  in  the  model.  The  size  of  the  plastic  zone,  w,  is 
determined  by  the  condition  that  the  stresses  should  be  bounded  at  the 
outer  edges  of  the  plastic  zone.  This  means  that  the  stress  intensity 
factors  due  to  the  external  loadings  (with  a replaced  by  a + w)  and  due 
to  the  ^ loadings,  should  sum  to  zero,  and  thus  for  the  two  problems 

considered  w is  the  solution  of: 


K (a  4-  w)  + K (a  + w)  = 0 

e p 


(4.13) 


-«r 


\ 


in  order  to  anihilate  the  stress  singularity.  The  Westergaard  Stress 
Function  for  the  problem  is 

Z(z)  = Zf (z-w,  a+w)  + Zp(z,w,a)  . (4.14) 


The  function  n c • i 1 1 in  general  oe  known  from  the  appropriate  elastic 
e 

solution.  ,D  may  oe  obtained  b\  applying  the  solution  for  wedging  force 
load?  or  a tract  as  a Green's  Function  to  generate  r.nt  solution  for 
distribute  loads  of  intensity  v>> . or  example,  irwin  [7]  has  provided 

the  ippropriate  westergaard  Stress  Function  (figui e -*.6  • which  is  usee 
as  3 Green's  Function  in  generating  L . Irwin  gives- 


Figure  4.6.  Wedging  force  on  a crack. 


By  replacing  a1  with  a + w.  shifting  the  origin  to  the  right  "a, 

and  replacing  P b\  - d:  with  the  distance  along  x from  the  new 
y s 

or i in  to  - , the  Green's  Function  for  generating  Z is  obtained 

vs  ? 

After  integrating  from  - = 0 to  : - v one  has 
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which  is  an  exact  solution  to  the  strip  yield  model,  since  no  assumption 
has  been  made  regarding  the  plastic  zone  size. 


By  assuming  that  the  yielding  is  small  scale,  a general 
be  obtained  for  the  plastic  zone  sizp  v in  terms  of  K and  ~ 

Sana t tor  (+.19)  nv  replacing  the  secant  tern  with  the  first 
of  its  expansion,  i e., 
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Assuming  w « a and  neglecting  all  w/a  terms  in  comparison  to  unity, 

one  has  K (a  + w)  ' K (a)  = K and  Equation  ('4.21')  reduces  to 
e e e 


w 
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(4.22) 


Recall  that  for  Mode  III,  R^  = ~ ■ r--1  . Thus,  again  one  sees  that  for 

t 

yp 

small-scale  plasticity  (i.e.,  small  zone  size),  the  plastic  zone  size 
depends  on  the  loading  and  geometry  only  through  the  elastic  stress 
intensity  factor. 

Equation  (4.22)  was  derived  from  stress  fields  for  cracks  in 
infinite  planes  These  results,  howeier,  will  be  valid  for  finite 
planes,  provided  that  the  plastic  zone  size  is  negligible  not  only 
in  comparison  to  crack  length  but  also  to  alt  planar  dimensions  of  the 
cracked  body  An  important  feature  of  the  strip  yield  model,  as  just 
illustrated,  is  the  fact  that  the  methods  of  arialvsis  are  essentially 
the  same  as  those  used  in  the  conventional  elastic  treatment  of  crack 
problems . 

4 u Limitations  of  the  Crack-Tip  Stress  Tield  Theory 

A certain  amount  of  nonlinear  deformation,  including  plastic 
yielding,  can  occur  without  grossly  affecting  the  surrounding  elastic 
crack-tip  stress  field.  Moreover,  such  nonlinearities,  if  embedded  in 
identical  stress  fields,  will  likewise  be  identical  and  self-compensating 
in  fracture  toughness  testing.  Thus,  it  is  important  to  determine  the 


e 


\ ..  - Vs  '*  - a It 


' v W 


As  r becomes  larger.  Equation  (4.23)  becomes  more  and  more  approximate. 
The  exact  stress  may  be  obtained  from  the  Westergaard  Stress  Function 
[see  Equation  (2.32)) 
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A comparison  of  Equations  (4.23)  and  f4 ,2~)  is  given  in  Figure  4 7 
xhe  increasing  inaccuracy  in  the  approximate  stresses  is  apparent  as  r 
increases.  It  will  be  useful  in  the  following  to  have  the  ratio  of  the 
approximate  and  exact  stresses.  Namely, 
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(4  25.1 


Consider  next  the  problem  shown  in  Figure  4.8.  Loadings  in  this 
and  the  preceding  problem  represent  opposite  extremes.  Westergaard 
has  given  the  following  stress  function  for  this  problem  (or  take  b = 0 
in  Equation  (4.15)), 
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From  Equation  (2  29),  the  stress  intensity  factor  for  this  case  is 
calculated  to  be 


(4.27) 


Substituting  into  Equation  (2.11),  the  approximate  stresses  are 
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I'ving  the  Westergaard  Stress  Function,  Equation  (4.26''  into  Equations 
(2.6),  the  exact  stresses  are  determined  to  be 
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These  two  equations,  (4.28)  and  (4.29),  are  compared  in  Figure  4.8. 


70 


JWKfc.yflhs*. t&SWfc  '/Mvm  .•*  % JL &?■.  \\  v*  - v—  ■ »•  * 


Figure  4.8  Comparison  of  approximate  and  exact  solutions 
for  wedging  forces  on  a finite  crack. 

Both  configurations  (Figures  4.7  and  4.8)  illustrate  the  accuracy 
of  the  stresses  when  i is  small  and  the  departure  from  accuracy  as  r 
increases.  For  the  two  problems,  results  indicate  similar  deviations 
of  the  approximate  stresses  from  the  exact.  Thus,  if  the  tolerable 

plastic  "one  si~e  can  be  determined  for  one  problem  it  sho.  1 ’ he  ippl 
cable  to  others 

While  for  brittle  problems,  the  usefulness  of  the  approximate 
stresses  for  predicting  fracture  behavior  is  not  questioned.  Figures 
4.7  and  4.8  contain  strong  implications  pertaining  to  certain  experi- 
mental techniques.  For  example,  in  using  photoelasticity,  moire,  etc. 
measurements  must  be  made  very  near  the  crack  tip  before  favorable 
comparison  with  the  approximate  stresses  could  be  expected. 

Returning  to  the  question  of  the  tolerable  plastic  zone  size, 
recall  the  results  of  the  Mode  III  solution  (the  strip-yield  results 
were  comparable), 
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in  Mode  I notation,  ('sins,  the  stress  intensity  factor  ^ - 
for  the  first  configuration. 


3 = ( — 


(-31) 


in  accordance  with  the  Mode  III  solution,  imagine  the  crack  tip  to  be 
located  at  the  center  of  the  plastic,  zone.  The  edge  of  the  plastic 
zone  will  be  located  at 


(4.32) 


Thus  the  size  of  the  plastic  zone  is  determined  by  the  ratio  a/a 

r ys 

The  linear  fracture  mechanics  approach  is  usually  assumed  to  be  valid 
up  to  stress  levels  of  0.8  of  the  yield  strength.  In  fact,  the  limi- 
tation placed  on  the  above  ratio  in  the  fifth  ASTM  Special  Committee 
Report  was 


< 0.8 


Substituting  the  upper  limit  of  stress  a = 0.8  a into  Equation  (4.32), 
the  relative  size  R / a for  reasonable  accuracy  is  about  0.3.  From 

y 

Equation  (4.25),  the  actual  stresses  deviate  about  20  percent  from  the 
approximate  stresses.  It  thus  appears  that  the  plastic  zone  at  a crack 
tip  may  be  fairl-  large  .••ithout  grossly  disturbing  tru  user.  lriss  M 
the  elastic  stress  field  approach  It  is  pointed  out  that  this  is  only 
a rough  estimate  intended  to  illustrate  the  usefulness  of  the  elastic 
stress  approach  in  the  presence  of  a given  amount  of  yielding  Require- 
ments on  the  plastic  zone  size  for  obtaining  valid  plane-strain  fracture 
toughness  data  are  much  more  strigent  than  the  above  would  seem  to 
indicate.  This  is  necessary  in  order  to  obtain  fracture  toughness 
values  which  are  independent  of  specimen  dimensions  (see  Chapter  5). 


Chapter  5.  FRACTURE  TOUGHNESS  TESTING 


The  purpose  of  fracture  toughness  testing  is  to  determine  the 
critical  stress  intensity  factor  for  a given  material  in  a given  environ- 
ment The  basic  procedure  is  to  load  a precracked  specimen  to  failure, 
observe  the  failure  load  and  crack  length,  and  use  the  appropriate  solu- 
tion for  the  rest  specimen  to  calculate  the  value  of  the  critical  stress 
intensity  factor.  Other  tests  are  aimed  at  determining  the  rate  of 
crack  growth  under  cyclic  loads  ar.d  the  fracture  toughness  under  long- 
term sustained  loads  in  a corrosive  environment.  Though  conceptually 
simple,  there  are  numerous  considerations  tf  basic  material  behavior  and 
experimental  procedure  which  must  be  accounted  for  ir.  order  to  obtain 
useful  fracture  data.  A comprehensive  collection  of  fracture  data  have 
been  tabulated  in  the  Damage  Tolerant  Design  Handbook  150’.  ASTM 
specifications  151]  now  control  the  method  of  test  for  plain-strain 
fracture  toughness  testing.  \ review  of  developments  leading  to  the 
adoption  of  the  ASTM  specitications  is  included  m the  ASTM  Special  itcr.ni- 
cal  publications  STP  410  152]  and  STP  463  [53].  Fatigue  crack  growth 
data  and  additional  fracture  toughness  testing  are  included  in  the 
publications,  STP  513  [54]  and  STP  514  [55],  Part  I and  Part  II,  respec- 
tively, of  the  proceedings  of  the  1971  National  Symposium  on  Fracture 
Mechanics  Since  test  methods  are  described  in  detail  elsewhere,  only 
a summarv  of  basic  ideas  is  presented  here. 

5.1  Characteristics  of  tne  Critical  Stress  Intensity  Factor 


because  of  its  technological  importance,  most  fracture  toughness 
testing  has  been  directed  to  the  Mode  T case  The  value  of  the  stress 
intensity  factor  at  which  crack  propagation  occurs  is  known  as  K . 

Irwin  [56]  determined  chat  (or,  alternatively,  the  critical  strain 

energy  release  rate  gc^  varied  with  plate  thickness  His  results, 

taken  from  Reference  [57],  for  tests  on  plates  of  7075-T6  aluminum  are 
shown  in  Figure  5.1.  The  curve  shows  a decreasing  toughness  for 

plates  or  increasing  thickness , while  the  fracture  surface  makes  a 
transi.ion  irom  full  slant  to  nearly  square  (Figure  5.2’  For  thin 
specimen-  the  stace  of  stress  at  tne  crack  tip  approaches  that  of  plane 

stress  (o  =0).  In  the  absence  of  the  state  of  stress  is  biaxial 

zz  zz 

permitting  the  formation  of  a relatively  large  plastic  zone  (see  Chapter 
4t.  The  ability  to  yield  thus  increases  the  fracture  toughness  Since, 
for  thin  specimens  the  state  of  stress  fs  approximately  biaxial, 

, the  maximum  shear  stress  occurs  on  a plane  at  approximately 

xx  yy  K 

-+5  degrees  to  the  cracK  plane  and  a slant  tor  V-slant)  fracture  results 
(i.e.,  P = 0)  as  illustrated  by  Figure  5.2  and  as  shown  by  the  dotted 
line  in  Figure  5.1.  Focusing  attention  on  the  crack  tip  for  this  case, 
it  can  be  seen  that  failure  occurs  by  a combination  of  opening  mode  and 
edge  tearing  or  a combination  of  Modes  I and  III.  As  the  specimen 


73 


..very*  ::#v* 


100 


1000 


P = 100X/B 


THICKNESS.  B (in.) 


gure  5.1.  Appearance  of  fracture  surface  and  effe 
of  place  tnickness  on  fracture  toughness 
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Figure  5.2  Types  of  fracture 
taken  normal  to  direction  of 


surfaces  il^ustratec 
crack  prooagation 


1 


thickness  increases,  the  state  of  stress  along  the  crack  front  approaches 

that  of  plain  strain,  o = ( + o ),  resulting  in  a state  of  tri- 

r zz  xx  yy  . 

axial  tension;  plastic  flow  is  inhibited  and  the  fracture  toughness 

decreases  while  the  fracture  surface  becomes  more  nearly  square.  A 

state  of  plane  stress  still  exists  on  the  plate  surface,  and  shear  lips 

form  there  as  shown  in  the  section  view  of  Figure  5.1.  Figure  5.2  gives 

a pictorial  representation  of  the  thin-to-thick  fracture  surface 

trans ition . 

As  the  specimen  thickness  is  increased  further,  the  curve  levels 

out  and  approaches  a constant  value.  This  value  is  designated  as  the 
plane-strain  fracture  toughness  (historically  , when 

determined  in  accordance  with  ASTM  specifications,  is  taxen  to  be  a 
material  constant  characteristic  only  of  the  given  materia) 

Following  Wilhelm  [45],  the  variation  of  the  critical  stress 
intensity  factor  with  thickness  for  several  materials  is  shown  in 
Figure  5.3.  The  values  are  nominal  values  only  and  are  shown  to 
illustrate  behavior.  Some  materials,  7075-T6  and  Ti-B120VCA,  quickly 
approach  a constant  value  of  while  others  require  greater  thick- 
nesses . 
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Figure  5.3.  Nominal  critical  stress  intensity  factors 
for  several  materials  as  a function  of  thickness. 
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5.2  Plane  Strain  Fracture  '''o^aness  Testing 


Numerous  specimen  designs  have  been  used  in  the  past  to  obtain 
fracture  toughness  data,  but  only  two  specimens  are  now  recognized  as 
ASTM  standards  [511  Thev  are  the  three-point -bend  specimen  and  the 
compact -tens ion  specimen.  These  are  shown  in  Figure  5.4.  Each  specimen 
contains  a crack-s carter  notch.  The  che  /eron  notch  is  shown  in  Figure 
5. 4,  but  several  forms  ar,  acceptable  The  notch  must  terminate  witi 
a fatigue  crack  of  at  least  0.05  inch  in  length.  The  fatigue  cracking 
Load  must  be  limited  to  a relative l-  low  value  to  insure  a suitably 
snarp  crack  with  a small  crack-tip  plastic  zone.  The  curvature  of  the 
fatigue  crack  front  must  also  be  within  certain  limits.  The  original 
crack  length  is  measured  after  fracture  at  three  locations-  fa;  at  the 
center  of  tne  crack  front,  (b  ■ midway  between  the  center,  and  < c)  the 
end  of  the  crack  front  on  each  side.  The  average  of  these  three  is  then 
used  as  tne  crack  length,  a,  ir,  subsequent  calculations  The  dimensions 
or  tne  specimens  are  controlled  by  the  crack-tip  plastic  zone  size.  The 


crack  length,  a,  and  the  thickness,  B,  must  exceed  2.5 


K_\2 


ys 


)■ 


where 


ys 


is  the  0.2  percent  offset  yield  strength  of  the  material.  This  is 


approximately  50  times  the  plane-strain  plastic  zone  size. 


The  specimen  is  loaded  to  failure  at  a specified  load  rate  either 
in  three-point  bending  or  tension  depending  on  the  specimen  configura- 
tion used.  During  loading,  a displacement  gage  which  spans  the  notch 
measures  the  crack  opening  displacement  (COD).  The  load  and  COD  are 
recorded  autograph ical ly,  resulting  in  a curve  similar  to  one  of  the 


tree  shown  in  Figure  ;.5.  The  load  used  to  calculate  a conditions 


value  v,  of  the  critical  stress  intensity  factor  is  called  P A 
0 * x 


secant  line  0P_  is  constructed  on  the  test  record  having  a siope  of 


95  percent  of  the  slope  of  the  tangent  OA  to  the  initial  linear  part  of 


the  record.  For  the  respective  type  of  record,  P^  then  becomes  the 


■■aiti  indicated  in  i igurt 


"a  tid 


test,  tne  rat  it 


P F,  must  not  exceed  I 10. 
max  Q 


..  conditional  critical  stress  intensity  tactor  1L  in  psi-in. 
is  calculated  as  follows; 


1 _ 


Three-point  bend  specimen  (Figure  5.4a): 
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Figure  5.5  Principal  types  of  load-displacement  records 


where 


load  from  load  - deformation  curve,  pounds 


B = specimen  thickness,  inches 


width  ( depth ) of  specimen,  inches 


b!  Compact  specimen  (Figure  5.4b) 


such  as  minimum  specimen  size  (in  comparison  with  2.5 


ment 


fatigue  cracking  load  level,  and  crack-tip  shape  does  the  test  become 
a valid  test.  Then  the  conditional  value  becomes  K , the  plane 


strain  fracture  toughness.  Plane-strain  fracture  toughness  testing  is 
unusual  in  that  the  test  requirements  (specimen  size,  etc.)  depend  upon 
K^c,  the  quantity  to  be  determined.  Therefore,  there  is  no  advance 

assurance  that  a valid  will  be  determined  in  a particular  test. 


- 


The  range  of  application  of  the  ASTM  test  method  is  restrictive.  The 
method  does  not  apply,  for  example,  to  brittle  materials  having  a thick- 
ness less  than  0.25  inch.  Also,  since  specimen  size  increases  proportional 


then  for  low-strength,  high-toughness  materials  the  method 


leads  to  specimens  of  a prohibitively  large  size.  Nor  is  it  possible 
to  determine  if  any  dimension  of  the  available  stock  of  material 


is  insufficient  to  provide  a specimen  of  the  required  size.  By  giving 
due  consideration  to  the  principles  involved  in  the  development  of  the 
standard  method,  it  is  sometimes  possible  for  the  investigator  to  design 
a useful  fracture  toughness  specimen.  For  a good  example  of  this  see 
Kendall  and  Hussain  [58], 


The  value  K^  determined  in  accordance  with  Reference  51  can  be 

considered  to  be  a material  constant.  It  can  serve  the  following 
purposes : 

a)  In  research  and  development  to  assess  the  effects  of  metallur- 
gical variables  such  as  heat  treatment  or  of  fabrication  operations 
such  as  welding  on  the  fracture  toughness  of  new  or  existing  materials. 

b)  In  design  applications  or  service,  evaluation  to  establish 
the  suitability  of  a material  for  a given  specific  application  where 
the  load  level  and  maximum  flaw  size  are  known. 

c)  For  specifications  of  acceptance  and  manufacturing  quality 
control . 


5 . 3 Plane-Stress  and  Transitional  Fracture  - Toughness  Testing 
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Currently,  no  ASTM  standard  exists  for  determining  plane-stress 
or  transitional  values  of  Kc-  Development  of  a standard  method  has 

been  hampered  by  the  degree  to  which  crack  tip  plasticity  and  associ- 
ated stable  crack  extension  occurs  prior  to  fracture.  Consequently, 
considerable  emphasis  must  be  directed  to  monitoring  the  slow,  stable 
tearing  portion  of  the  fracture  process  to  determine  an  effective 
critical  crack  length.  A number  of  experimental  guidelines  have  been 
developed,  as  explained  in  Reference  50,  and  were  used  as  a criteria 
in  selecting  K^  data  for  inclusion  in  that  handbook. 


t 
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In  1970  a research  program  was  initiated  at  the  Naval  Research 
Laboratories  (NRL)  for  developing  a standard  testing  procedure  for 
determining  K values  for  thin  sheets  [59],  This  research  has  resulted 


in  a proposed  standardized  test  for  screening  purposes  using  the  center- 
cracked  tensile  (CCT)  sheet  specimen.  Some  of  the  details  of  this 
method  differ  from  those  summarized  in  Reference  50  Because  of  the 
extensive  research  conducted  m arriving  at  the  method  proposed  in 
Reference  [59],  the  general  features  of  that  method  are  summarized  here 

The  CCT  specimen  of  thickness  B and  width  W is  shown  in  Figure  5.6. 
The  recommended  length  L is: 

L = 2W,  W <12  in. 

L = 1.5  W,  W > 12  in 

where  L is  measured  between  the  innermost  grip  pins  (grips  are  detailed 
in  Reference  [59]).  The  recommended  final  crack  length  2ac  is  approxi- 
mately 0.333  W.  (The  original  crack  length  2a°  will  be  smaller.)  The 

starter  slot  ideally  should  be  subtended  by  fatigue  cracks  although  in 
some  cases  a slot  with  a root  radius  of  approximately  0.001  to  0.002 
inch  will  effectively  simulate  a fatigue  crack.  For  specimens  thinner 
than  0.123  inch,  crack  buckling  is  suppressed  by  the  use  of  rigid, 
lubricated  face  plates  attached  to  the  central  portion  of  the  specimen. 


To  determine  the  dimensions  of  the  specimen  the  plastic  zone  size 
is  estimated  as 
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During  the  fracture  test  a COD  gage  is  used  to  measure  crack 
opening  (similar  to  plane-strain  testing)  as  the  load  is  applied.  The 
effective  crack  length  2a  is  determined  after  the  test,  not  by  mea- 
suring but  from  a compliance  calibration  curve.  To  determine  the  com- 
pliance calibration  curve,  specimens  are  prepared  with  slot  lengths 
varying  through  a range  of  values  of  2ao^w  such  as  0.1,  0.2,  0.3,  0.4, 

0.5,  0.6  The  specimens  are  loaded  elastically  while  P and  COD  are 
graphed  by  an  x - y recorder  (Figure  2.1).  This  results  in  a compliance 
calibration  curve  which  can  be  plotted  nondimensionally  as  shown  in 
Figure  5.7  By  measuring  the  quantity  EBtC0D]/P,  where  E is  Young's 
modulus,  during  the  actual  fracture  test  one  can  then  use  the  compli- 
ance curve  (Figure  5 7)  to  determine  the  effective  critical  crack 
length  a . This  procedure  eliminates  the  problem  of  physically  measuring 

the  crack  length  at  the  point  of  instability.  The  need  for  adding  a 
piastic  zone  correction  to  the  crack  length  is  also  bypassed  since  the 
conpliance  curve  gives  the  effective  critical  crack  length. 


2a/ W 


Figure  5.7.  Example  compliance  calibration 
curve  for  a center  - cracked  tensile 
specimen. 

Figure  5.8  shows  two  types  of  load-displacement  curves  for  a plane- 
stress  fracture  test.  In  curve  A,  the  curve  rises  linearly  to  a given 
point  before  crack  growth  starts,  then  continues  to  rise  with  accompany- 
ing crack  growth  until  final  instability  occurs.  After  the  maximum 
load  is  reached  on  Curve  B,  the  crack  continues  to  grow  under  constant 
load  until  final  separation. 
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Figure  5.8.  Load  - COD  curves  showing  two  types 
of  crack  growth  response. 

To  find  Kc , both  curves  are  analyzed  in  the  same  way.  A straight 

line  (a)  is  drawn  through  the  linear  portion  of  the  curve.  The  value  of 
EB[COD] /P  is  calculated  and  compared  with  the  value  known  for  2aQ/w  from 

the  compliance  curve  (Figure  5.7),  and  the  difference  is  noted.  There 
will  be  a slight  difference  due  to  variations  of  alignment,  etc.  A 
secant  line  (b)  is  then  drawn  to  the  first  point  of  maximum  load  value; 
EB[COD]/P  is  calculated  for  this  line  and  corrected  by  the  difference 
noted  above.  The  corrected  value  of  EB[COD]/P  is  then  used  on  the  appro- 
priate calibration  curve  (Figure  5.7)  to  determine  the  effective  critical 
crack  length  Finally,  «c  is  calculated  from  the  following  formula: 


r 

/2a  \ 

/2arV 

K = o /a~ (1.77) 
c ' c 

1 - »■>(-)  < 

(nr)  j 

I'ne  aDove  is  only  a summary  of  plane  stress  testing,  and  for  furthei 
details  the  reader  is  referred  to  the  handbook  [50],  and  the  paper  by 
Sullivan  and  Stoop  [59]. 

5 . U Fatigue  Crack  Propagation 

The  rate  of  growth  of  a fatigue  cracs  under  a time-varying  load 
is  influenced  by  a number  of  variables,  including  stress  level,  crack 
length,  specimen  geometry,  load  frequency,  material,  environment,  and 
temperature.  The  effects  of  stress  level,  crack  lengtns,  and  specimen 
geometry  can  be  properly  accounted  for  by  means  of  the  stress  intensity 
factor.  For  a fatigue  loading,  the  applied  stress  varies  with  time  and 
since  the  stress  intensity  factor  is  proportional  to  stress  then  the 
stress  Intensity  factor  exhibits  the  same  kind  of  variation  with  time. 
That  is,  K has  the  form 
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(5.8) 


For  a sinusoidal  load  the  crack  extends  an  amount  La  during  each 

load  cycle  _.N,  where  N is  the  number  of  cycles.  In  fatigue  testing  the 

incremental  crack  growth  rate  ^a  . N is  determined  as  a function  of  the 

stress  intensity  factor  range  where  ,K  = K - K , and  the  other 

° max  min 

variable^,  such  as  temperature,  are  normally  treated  as  parameters, 
constant  for  a given  set  of  data  The  quantity  a A_At,  determined  by- 
monitoring  crack  length  and  the  number  of  cycles,  is  generally  written 
as  a continuous  quantity'  da/dN.  The  data  are  presented  as  a log-log 
plot  of  da/dN  (or  d(2a)  dN)  versus  ,K  as  shown  by  Figures  5.9  and  5. 1C, 
taken  from  Reference  50 

For  a wide  range  of  IK  values  the  data  approximate  a straight  line 
Noting  this,  Paris  suggestec  a power  law  of  the  form  (for  example,  see 
References  60  and  61), 


di = c(JiK)n 


where  c depends  upon  the  load  ratio,  R 


(5.9) 


c . /c  , 

min  max 


the  material,  and 


other  test  parameters.  Instead  of  the  stress  ratio,  Paris  used  the 
relative  mean  load  r defined  in  the  present  notation  by  y - 


j /A a . 

mean 


! 

§ - 

> 

where 

I 

-1? 

C * 1.0  x 10 

■ 

; 

C = 1.6  x 10‘12 

1 

/ 

1 

C - 2.0  ■ 10"12 

For  several  materials  with  y = 0.5  (R  = 0),  a value  of  n = 4 gave  a 
reasonably  accurate  fit  to  the  data,  as  shown  by  Figures  5.11  (a),  (b ; , 
and  (c)  from  Reference  60.  In  fact,  based  on  Figure  5.12  for  steels 
loaded  with  7 - 0.5,  Irwin  and  Paris  [ 62 1 , recommended  for  estimation 
purposes 


(:  .10) 


steel ) 


and 


all  steels') 

da/dN  = inches  per  cycle 
K = ksi-inch^2  . 
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STRESS  INTENSITY  FACTOR  RANGE,  (k«i-in.1/2  ) 

Figure  5.9.  Variation  of  fatigue-crack-growth  rate 
with  K for  2024-T3  alloy. 
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STRESS  INTENSITY  FACTOR  RANGE.  AK  (ksi— in. 1 /2 > 

Figure  5.10.  Variation  of  fatigue-crack-growth  rate 
with  K for  7075-T6  aLLov 
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:a'  DATA  TOR  7075- T6  ALUMINUM  ALLOY 


Figure  5.11.  The  exponential  equation  of  Paris 
fitted  to  the  data  for  various  materials 
(Figures  from  Reference  60). 


87 


* 


L u.  I 


Figure  5.12.  Crack-growth  rate  da/dN  as  a function 
of  fatigue  cycle  range  iaK.  A302B  data  (from 
Reference  62 ) 

Equation  (5  10')  is  a crude  relationship,  not  containing  the  effects 
of  load  ratio,  frequency,  environment,  or  temperature,  etc.  Hence, 
it  is  only  recommended  for  estimates  of  whether  fatigue-crack  growth 
is  a real  problem  or  not  for  an  anticipated  flaw,  and  not  recommended 
for  actual  estimates  of  tne  fatigue  life  of  structures. 

Recent  results  by  Barsom  [63]  indicate  that,  at  a low  stress 

intensity  factor  range,  the  rruck  propaga: icu  rate  in  high  yield 

trength  steels  conform  more  closely  to  a second  power  law  than  to 
a fourth  power  law.  The  equation  recommended  by  Barsom,  which  was 
verified  for  19  types  of  steels,  is 

= 0.66  10'8  (_JC)2-2:i  (5  11 

dN  v 

, -9 

tor  tne  steel  D6A0 , values  oi  n = 2 35  and  C = 1 . 9 ’ 10  were  recent! 

f > 'nd  bv  Paris  [641  Much  earlier,  Liu  [65]  had,  in  fact,  claimed, 
based  on  deduction,  a value  of  n = 2.  It  seems  clear  though  that  the 
data  fail  to  support  any  fundamental  value  of  n.  Thus,  Equation  (5  9 
must  continue  to  be  viewed  as  an  empirical  relationship,  useful 
nevertheless . 
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f 

Equation  (5.9)  does  not  account  for  two  important  effects:  the 

variation  in  the  crack  growth  rate  due  to  the  load  ratio,  R,  and  the 

1 instability  of  the  crack  growth  when  the  stress  intensity  factor 

approaches  the  Kc  value  for  the  material.  To  account  for  these,  Forman 

et  al.  [66 J proposed  the  following  equation 

da  . C QuK)n  . ( 

dN  (1  - R)  K - AK  v ' 

For  2024-T3  and  7075-T6  aluminum  and  n = 3,  Forman's  equation  gave  good 
correlation  with  data  for  several  values  of  R.  For  the  same  two  materi- 
als, Hudson  [6"’l  further  verified  Forman's  equation  by  showing  that  it 
produced  an  excellent  fit  for  a wide  range  of  R ratios.  In  theory,  for 
any  problem  for  which  there  is  an  analytical  expression  for  K,  Equation 
(5  12)  can  be  integrated  to  determine  the  crack  length  (or,  alterna- 
tively, K)  after  a given  number  of  cycles.  In  practice,  for  most  cases, 
numerical  integration  must  be  resorted  to  due  to  complex  expressions 
for  K,  especially  if  the  fatigue  load  is  nonuniform,  i.e.,  if  R and  Ao 
are  not  constant. 


As  a special  case,  Forman  illustrated  the  case  of  an  infinite  piate 
containing  a crack  of  length  2a.  Then 

AK  = Aa  y^ra  . (5.13) 
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While  equations  such  as  (5.9)  and  (5.12)  represent  trend  of 
fatigue  data  it  must  be  kept  in  mind  that  environmental,  material, 
temperature,  and  frequency  effects  all  are  lumped  into  the  constants 
C and  n.  Moreover,  for  some  materials  and  loading  conditions,  other 
phenomona  require  attention.  Paris  [68]  has  investigated  the  idea  of 
threshold  stress  intensity,  value  of  AK  below  which  the  crack  propaga- 
tion rate  becomes  vanishingly  small.  For  values  of  R <.  0,  crack 
closure  can,  for  some  materials,  cause  additional  fatigue  damage 
resulting  in  increased  crack  growth  rates.  This  has  been  observed 
by  Elber  [69]  and  Hudson  [67]  among  others.  For  additional  information 
on  recent  fatigue  crack  testing,  the  reader  is  referred  to  ASTM  STP 
486  [70],  and  ASTM  STP  513  [54].  For  those  interested  specifically  in 
D6AC  steel  the  report,  "Crack  Behavior  in  D6AC  Steel  - An  Evaluation 
of  Fracture  Mechanics  Data  for  the  F-lll  Aircraft"  [71],  contains  the 
largest  amount  of  crack  growth  data  ever  generated  for  any  one  material 
application . 

For  estimating  flaw  life  with  accuracy  one  can  select  appropriate 
data  from  the  Damage  Tolerant  Design  Handbook  [50]  or  resort  to  labora- 
tory tests  to  measure  da/dN  versus  AK,  under  the  expected  conditions  of 
load  ratio,  environment,  frequency,  and  temperature.  A log-log  plot  of 
da/dN  versus  AK  then  gives  the  basic  crack  behavior  of  the  material. 
Numerical  integrat ion  (See  Reference  [ 66 ])  or  the  report  by  Engle  [72]) 
of  the  crack  growth  rate  from  the  initial  flaw  size  to  the  critical 
flaw  size  is  then  required  to  estimate  the  fatigue  life  of  the 
structure . 

5.5  Stress  Corrosion  Cracking 

Under  static  load  and  the  influence  of  an  aggressrre  or  corrosive 
environment,  stable  crack  growth  can  occur  at  values  of  K^.  considerably 

below  the  critical  value.  Such  slow  crack  growth  is  referred  to  as  sub- 
critical  crack  growth.  It  can  eventually  lead  to  a critical  flaw  size 
and  failure.  Generally,  two  types  of  tests  are  conducted  to  assess  the 

eii'e ct  o : an  aggressive  environment,  (a)  cime-to-failure  tests,  and 
(b)  sustained  crack  growth  rate  tests.  Numerous  types  of  specimens  are 
used,  including  those  already  covered  in  this  chapter.  Details  of  speci- 
mens as  wel]  as  data  for  several  materials  are  included  in  Reference  50. 


In  time-to-failure  tests,  precracked  specimens  immersed  in  the 
corrosive  environment  are  loaded  at  various  levels  of  stress  intensity 
and  the  time  to  failure  is  recorded  The  stress  intensity  is  reduced 
in  steps  until  a threshold  level  of  K^.  is  reached,  below  which  no 

failure  occurs.  This  K level  is  called  K . Values  are  tabulated 

1.  .Lo\_<U 

in  the  Damage  Tolerant  Design  Handbook 


The  difference  in  the  environmentally  determined  threshold  stress 
intensity  ( ^ ancl  the  plane-strain  fracture  toughness  (K^)  in  air 
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provides  a measure  of  the  susceptibility  of  the  material  to  stress 
corrosion  cracking.  The  trend  of  data  for  such  tests  as  indicated  by 
Wilhelm  [45]  is  shown  in  Figure  5.13.  In  high  strength  materials, 

Kj  c may  be  small  compared  to  while  in  lower  strength,  tougher 

materials,  may  approach  or  equal 


i 


Some  K_j.scc  data  have  been  obtained  in  natural  or  simulated  natural 

environments.  Most  data,  though,  have  been  obtained  in  aqueous 

solutions  of  3.5  weight  percent  NaCl.  This  solution  contains  about  the 
same  percentage  of  salt  as  seawater.  Other  environments  include  dis- 
tilled water  and  air  at  or  near  100  percent  relative  humidity 


Environmentally  sustained  crack-growth  rates  (da./dt)  have  been 
obtained  with  a variety  of  specimens  In  general,  results  are  obtained 
in  most  of  the  same  environments  as  the  KISCC  data.  Although  other 


methods  are  used,  crack  extension  is  usually  determined  by  optical  mea- 
surements of  the  crack  length  as  the  crack  propagates.  An  example  set 
of  results  from  a crack-growth-rate  test  is  shown  in  Figure  5.14. 
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F igure  5.13. 
cracking  in 


Delayed  failure  due  to  stress  corrosion 
an  aggressive  environment. 


91 


j 


Chapter  6.  DESIGN  APPLICATIONS  OF  FRACTURE  MECHANICS 


In  design,  the  fail  safe  philosophy  requires  that  during  the 
service  life  of  a structure,  any  fatigue  cracks,  defects,  or  damage 
will  not  lead  to  failure  prior  to  detection  during  regular  inspection 
■trnds.  This  means  chat  a structure  containing  readily  detectable 
damage  or  cracks  will  retain  sufficient  residual  strength  to  provide 
safe  seivice  upon  subsequent  load  history.  Thus,  concern  is  for  the 
remaining  strength  of  a structural  member  in  the  presence  of  an  exist- 
ing crack.  Fracture  mechanics  provides  the  technique  for  making  such 
strength  caiculations . Coupled  with  nondestructive  testing  and  reliable 
inspection  techniques,  fracture  analysis  can  help  ensure  the  integrity 
of  an  engineering  structure. 


Several  example  design  problems  follow.  Some  are  concerned  with 
only  residual  strength  calculations; others  involve  fracture -re si  stance 
comparisons  and  selection  of  materials.  In  the  illustrations,  frequent 
use  will  be  made  of  the  crack  solutions  which  appear  in  Chapter  7. 

Some  of  the  problems  have  been  selected  from  Wilhelm  [45]  and  Wessel  [73] 
others  have  been  devised  under  the  guidance  of  these  publications. 

The  problems  are  for  illustrative  purposes  only,  and  should  be  used 
only  as  guides.  Each  engineering  problem  will  have  its  own  unique 
features  and  soLation.  The  values  of  fracture  toughness  used  here  may 
not  correspond  to  the  proper  values  required  for  a particular  problem 
and  environment.  One  should  select  the  appropriate  fracture  toughness 
data  such  as  and  da/dn  from  the  Damage  Tolerant  Design  Handbook 

[50]  or,  if  this  is  not  possible,  obtain  the  data  from  a testing  pro- 
gram as  outlined  in  Chapter  5. 


o.l  Residual  Strength  of  a Center-Cracked  Panel  [45] 


Ir  is  assumed  here  t <at  a simple  rectangular  aluminum  panel  uni- 
formly loaded  in  tension  across  the  width  (Figure  6.1)  will  develop  a 
central through-the-thickness  crack  4 inches  long  before  detection  by 
the  established  inspection  routine  '’it  is  noted  that  althourd  this 
crack  geometry  occurs  infrequently,  it  provides  a good  illustrative 
example).  The  material  is  7075-Tb  aluminum  having  a thickness,  B,  of 
0.100  inch.  The  problem  is  to  determine  if  the  panel  will  fail  at  an 
operating  stress  of  45  ksi  before  crack  detection.  It  is  decided  to 
analvse  tne  panel  using  a plane- stress  toughness  value.  Figure  5.3 
indicates  that  for  this  thickness  K„  ~ 85  ksi- v/ in. 


The  solution  for  the  panel  is 


(6.1) 
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From  Table  3.1,  f ~ 1.03  for  this  aspect  ratio.  Then 


83  ksi  v/in. 


' jta  f 

c 


(i) 


= 33  ksi 


n(2  in. ) 1.03 


This  critical  stress  would  be  reduced  by  20%  to  30%  if  crack  buckling 
taK.es  place.  Thus,  the  sheet  can  sustain  a gross  area  stress  which  is 
less  than  33  ksi  in  the  presence  of  a 4-inch-long  crack.  Because  this 
is  277„  below  the  required  operating  stress  (45  ksi) , the  requirements 
are  not  met. 

From  this  illustration,  a change  in  inspection  schedule  to  detect 
a crack  smaller  than  4 inches  or  a reduction  of  operating  stress  level 
or  both  is  indicated. 

6 . 2 Strength  of  a Bracket  with  Cracks  Emanating  from  a Hole 

A proposed  design  for  a bracket  contains  a 0.5-inch-diameter 
hole  (Figure  6.2).  It  is  anticipated  that  a fatigue  crack  will  grow 
from  the  hole  to  a length  of  0.25  inch  before  detection  by  inspection 
The  bracket  is  loaded  with  a remote  stress  of  52  ksi.  The  material 
data  follow: 

a)  Material:  4340,  oil  quenched,  tempered  at  500*F 


b)  Material  yield  strength 

c)  Material  fracture  toughness: 


vs 


= 238  ksi 


K = 46  ksi-in. 
1 o 


1/2 


Can  the  design  be  considered  to  be  a safe  one7  To  answer  this  question 
the  critical  stress  will  be  calculated.  It  will  be  conservative  to 
assume  that  a crack  develops  from  both  sides  of  the  hole.  The  problem 

then  becomes  the  one  shown  in  the  inset  of  Figure  6.2. 

Any  stress  intensity  factor  solution  for  this  problem  must 
reflect  both  the  presence  of  the  free  surface  of  the  hole  and  the  finite 
width  of  the  bracket.  For  cases  such  as  this,  an  approximate  solution 
can  often  be  developed  by  applying  separate  correction  factors  for  the 


various  free  surfaces.  Let  f 

'a. 


hole  and  f. 


ft) 


id) 


be  the  correction  factor  for  the 


be  a correction  factor  for  the  finite  width  of  the 


plate.  The  quantity  a is  an  effective  crack  length  for  computing  the 


Kj  in  the  form 


finite  width  correction.  Apply  corrections  for  both  effects  by  writing 
rm 

i 

^fi(i)f2(f)  * (6-2) 


KI  ^ 
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Figure  6.2.  Brackec  containing  a hole. 


This  is  clearly  an  approximation  since  it  neglects  any  interaction  of 
the  two  effects.  Such  a procedure  is  common  practice  and  is  normally 
sufficiently  accurate  for  design  calculations. 

a 

Returning  to  the  problem  at  hand,  compute  the  ratio  — 


3.  _ 3 • 25  inch 
R 0.25  inch 


From  Bowie's  solution  (Figure  7.17) 


For  computing  the  finite  width  correction,  an  effective  crack  length  ag 
of  a + R will  be  used.  Thus, 


e a + R 0.25  in.  + 0.25  in. 


b b 


From  Table  3.1, 


3 in. 


= 0.167 


m ■ 1 


02 


Solving  Equation  (6.2)  for  the  stress  and  letting  critical  conditions 
apply,  one  has 


■c  £l(f)  h(f) 


46  ksiV  in. 

/ ir( 0 . 25  in.)  (1.45)(1.02) 


35  ksi 


Since  the  critical  stress  is  well  below  the  imposed  load  of  52  ksi,  the 
design  will  not  tolerate  an  inspectable  crack. 

It  is  interesting  to  investigate  the  maximum  crack  length  which 
could  be  tolerated  by  this  design.  Solve  Equation  (6.2)  for  the  crack 
size  and  assume  critical  conditions  at  the  imposed  stress  of  52  ksi  . • 


a = — 

c n 


‘i(«>  hHP) 


This  givres  the  solution  for  the  maximum  crack  length.  However,  the  hole 
correction  f^  and  the  width  correction  f?  are  both  a function  of  crack 

length.  For  this  crack  geometry,  estimates  can  be  made  for  f and  f . 

a + R L L 

As  — g > 0,  f^  1.00.  Assume  f^  ~ 1.00.  For  the  hole  correction 

a 

factor,  f.  , assume  ^^3.39,  which  corresponds  to  — = 0.  After  calculat- 
11  R 

ing  a^  one  can  determine  revised  estimates  for  f^  and  f and  iterate: 


I 


M#*  - 


= I 

ir 


46  ksi  ./in. 


c it  (52  ksi) (3.39) (1.00) 


= 0.0217  in. 


check 


a 0.0217 


R 0.25 


= 0.088 


for  which,  from  Figure  7.17, 


f = 2.81 


a + R 0.0217  + 0.7.5 


• = 0.091 


for  which,  from  Table  3.1 


f2  = 1.00  . 


Calculating  again 


i 46  ksi  yin. 

ac  " Jt  (52  ksi)  (2 . 81)  (1.00) 


= 0.032  in. 


Three  more  such  iterations  finally  lead  to  the  following  critical  crack 
length:  ^ 


a = 0.038  in. 
c 


This  is  a very  short  crack  length  ana  could  perhaps  correspond  to  a 
nick  caused  by  a tool  at  the  hole  edge;  therefore,  a re-evaluation  of 
the  design  is  inevitable. 


In  a situation  similar  to  this,  Wilhelm  [45]  lists  four  possible 
design  modifications.  They  are: 


a)  Reduce  allowable  operating  stress. 

b)  Reduce  hole  diameter. 


c)  Use  a different  material  (higher  K ) 

J-  L. 


d)  Reduce  plate  thickness  and  still  meet  stress  criteria. 


WZ  S»r  .< 


i 
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Functional  constraints  may  prohibit  reducing  the  hole  diameter.  The 
fourth  alternative  would  be  possible  for  some  materials.  By  reducing 
the  plate  thickness,  one  takes  advantage  of  the  material's  ability  to 
deform  plastically  which  results  in  higher  values  of  fracture  toughness. 
For  the  present  material  though,  it  is  felt  that  a moderate  reduction 
in  thickness  would  not  significantly  improve  the  fracture  toughness, 
and  for  a thickness  reduction,  a wider  bracket  would  be  required  to 
meet  stress  criteria.  In  this  case,  probably,  the  best  solution  is 
to  select  a different  material  having  a higher  toughness,  which  also 
mav  have  a lower  operating  stress  due  to  a lower  yield  strength.  This 
attacks  the  real  weakness  of  the  present  design,  i.e.,  a low  toughness 

steel  (K  = 46  ksi  N/in. ) with  a high  yield  strength  (o  = 238  ksi) 

Ivj  y s 

which  permitted  a high  operating  stress. 

6 , 3 Cracked  Bar  Loaded  in  Bending 

An  arm,  loaded  in  bending  as  shown  in  Figure  6.3,  is  assumed  to 
develop  a crack  as  shown.  A force  P of  either  2000  or  4000  pounds  is 
applied  at  the  center.  Two  aluminum  alloys  are  being  considered: 

2219-T851  (K  = 35.8  hsi/HT.)  and  2124-T851  (K  = 28.8  ksi  yi^T)  . 

I v IC 

If  a 1/4-inch  deep  crack  can  be  detected  during  inspection,  which  of 

the  two  materials  would  be  the  most  fracture  resistant  for  the  two 

loads  under  consideration? 

The  stress-intensi  ty  factor  solution  for  this  problem  is  shown 
in  Figure  7.46  for  span  length-to-depth  ratios,  s/b,  of  4 and  8,  and 
also  for  the  case  of  pure  bending.  Also  note  that  this  solution  is 
for  a beam  of  unit  thickness.  The  solution  is  (allowing  for  the 
thickness  B) : 


where  M is  the  bending  moment.  In  this  case 


£ _ 12  in . 
b 2 in. 


a _ 0.25  in. 
b 2 in. 


0.125 
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Calculate  K^. 

K = .6  .(6.Q°0— in-_ 1.b.).-_  v/jt(0.25  in.)  (0.99) 
(0.75  in.)(2  in.) 

= 10.53  ksi  - v/in7  (for  P = 2000  lb) 

and 


K = 21.2  ksi  - ynr.  (for  P = 4000  lb) 

Compare  these  calculated  values  of  the  stress-intensity  factor  to  the 
plane-strain  fracture  toughness  values. 


Material 

K c (ksi  /in.) 

~ — (2000  lb) 
1C 

KI 

ha 

(4000  lb) 

2219-T851 

35.8 

0.29 

0.59 

2124-T851 

28.8 

0.36 

0.73 

For  the  2000-pound  load,  either  of  the  two  materials  might  be  considered 
acceptable  since  the  K /K  ratio  is  fairly  low  for  both.  For  the 

i.  1L 

4000-pound  load,  a greater  difference  is  apparent  and  the  2219-T851 
alloy  would  be  preferred. 

6.4  Evaluation  of  Materials  Based  on  K Data  and  Design  Requirements  - 
Surface  Flaw  [73] 


Three  materials  — a steel,  aluminum,  and  titanium  - are  to  be 
compared  for  a particular  application.  One  of  the  primary  considerations 
is  minimum  weight,  and  the  yield  strengths  of  the  three  materials  have 
been  chosen  so  that  their  strength-to-weight  ratios  are  equivalent. 

The  allowable  stress  is  to  be  one-half  of  the  yield  stress.  The 
material  properties  are  as  follows: 


Alloy 

Density 

(lb/in.3) 

0 

ys 

(ksi) 

o ^Density 
(x  1000  in. ) 

KIC 

(ksi  yin. ) 

Steel 

0.234 

250 

880 

100 

Aluminum 

0.098 

85 

870 

30 

Titanium 

0.163 

140 

860 

80 

101 


'.v.  •v*-  >•  fiv  a r.,- 


Figure  6.4.  Semielliptical  surface 
flaw  to  be  used  in  a materials 
evaluation. 


The  basic  solution  for  the  surface  flaw  is  the  Irwin  Solution 
Equation  (3.33).  This  solution  has  been  widely  used,  both  for  analysis 
and  for  fracture  toughness  testing.  It  does  not  contain  a correction 
for  the  back-side  free  surface  of  the  plate,  and  hence  is  not  valid 
if  the  flaw  is  "deep;"  i.e.,  if  the  flaw-depth- to- plate- thickness 
ratio  is  greater  than  0.5.  Consequently,  numerous  corrections  to 
Equation  (3.33)  have  been  developed,  some  based  on  results  of  tests, 
others  based  on  analytical  solutions.  The  Shah -K obayashi  solution  [74] 
(Figure  7.74)  is  an  analytical  solution  which  applies  to  a wide  range 

of  flaw  aspect  ratios  and  flaw-depth-to-plate- thickness  ratios. 

It  ssentially  agrees  with  the  Irwin  Solution  for  long,  shallow  flaws. 

It  should  be  noted  that  the  Irwin  solution  classically  contains  a 

2 

plasticity  correction  (the  0.212  — — term)  which  is  not  included  in 

a 

ys 

Figure  7.74.  For  precision,  one  should  use  the  Shah-Kobayashi  solution; 
however,  because  of  its  simplicity,  the  Irwin  solution  will  be  used  here 
for  illustrative  purposes.  The  Irwin  solution  is: 


1.1  a Jn  a 


E - 0.212 


' VS  ' 


2 I 1/2 


To  calculate  the  critical  flaw  size,  solve  the  above  for  a and  let 
critical  conditions  apply: 


I 


- 


The  following  results  are  tabulated  for  the  three  materials: 


Material 

a. . . (=  0.5a  ) 

All  v ys 

(ksi) 

KIC 

(ksi  ,/in. ) 

a 

cr 

(in.) 

Length 

(in.) 

Steel 

125 

100 

0.212 

1.06 

Aluminum 

42.5 

30 

0.165 

0.83 

Titanium 

70 

80 

0.432 

2.16 

For  che  conditions  imposed,  the  titanium  alloy  is  the  most  fracture 

resistant.  This  conclusion  could  have  been  reached  by  considering  the 

ratio  K .T„/a  for  the  following  materials: 

IC  ys 


Material 


Steel  0.400 
Aluminum  0.353 
Titanium  0.572 


0.160 

0.125 

0.327 


The  titanium,  having  the  highest  K /c 


ys 


ratio,  could  be  expected  to  be 


the  toughest  material.  Since  the  applied  stress  is  a fixed  percentage 

of  a for  all  the  materials,  the  critical  defect  size  will  be  propor- 
ys  ’2 

tional  to  the  ratio  (K^/a  s)  • AH  other  terms  in  Equation  (6.4)  are 

2 

constant,  therefore  (K^/a ) provides  a good  index  to  the  relative 


critical  defect  size. 


Assume  that  with  inspection  techniques  the  minimum  detectable  flaw 
is  0.15- inch  deep  by  0.75- inch  long.  Calculate  the  critical  stress  for 
fracture  from 


/1.21  Jta 


This  is  conveniently  done  by  iteration  (because  of  the  0.212  a /a  g 
term),  which  converges  very  rapidly.  The  results  are: 
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c 

(ksi) 


Material 

Steel 

Aluminum 

Titanium 


Note:  These  values  differ  slightly  from  those  in  Reference  73. 

With  this  information,  the  designer  can  now  re-evaluate  his  original 
choice  of  a design  stress  equivalent  to  one-half  of  the  yield  strength. 
With  assurance  that  a defect  of  the  type  used  in  the  foregoing  can  be 
reliably  detected  in  the  component,  the  designer  can  proceed  to  estab- 
lish the  safety  factor  in  his  design.  The  steel  looks  better  at  first 
since  it  requires  the  highest  stress  for  fracture  in  the  presence  of 
a flaw.  To  assess  the  safety  factor,  however,  one  must  consider  the 
ratio  of  the  critical  stress  to  the  design  stress  (1/2  a ).  Results 
are:  ys 


a 

a 

a 

ys 

ALL 

c 

c 

Material 

(ksi) 

(ksi) 

(ksi) 

°ALL 

Steel 

250 

125 

148 

1.18 

Aluminum 

85 

42.5 

45 

1.06 

Ti tanium 

140 

70 

115 

1.64 

Thus,  the  titanium  provides  the 

greatest 

safety 

factor  and  resistance 

to  fracture.  Steel  and 

aluminum 

are  marginal. 

The  designer  could 

alter  the  design  stress 

for  each 

of  the  materials  to  provide  whatever 

safety  factor  is  desired 

. Any  decrease 

in  the 

design  stress  to  raise 

the  factor  of  safety,  however,  would  result  in  an  increase  in  weight. 

6.5  Eccentric  Crack  Growing  from  a Rivet  Hole  [45] 

A crack,  starting  at  the  rivet  hole  (Figure  6.5)  has  grown  to  a 
length  of  3 inches  prior  to  being  detected.  The  panel  is  made  of 
7075-T6  aluminum  ^ith a plane-stress  fracture  toughness  of  approximately 

85  ksi  - v/in.  The  operating  stress  is  30  ksi.  What  are  the  critical 
conditions  for  fracture  for  this  crack  geometry  and  loading  and  should 
the  panel  be  replaced7 
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Figure  6.5.  Crack  growing  from  an  edge 
rivet  hole. 


In  solving  this  problem,  the  restraining  effect  of  the  adjoining 
panels  or  stiffeners  to  which  the  panel  is  riveted  will  be  ignored. 
(Panels  with  stiffeners  will  be  considered  in  the  next  example.)  This 
problem  can  be  treated  in  an  appropriate  manner  bv  simplifying  it  to  a 
panel  containing  an  eccentric  crack  of  length 

2a  ~ a + d = 3.1875  in. 
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for  a 3 16-inch  rivet  hole.  Figure  7.39  provides  the  solution  for  this 
problem  which  includes  both  the  eccentricity  and  finite  plate  width 
effect, 


4 = a 


Jta  F 


(v  Vi 


(6.5) 


where  b is  plate  half-width  and  e is  eccentricity.  A correction  factor 
F is  provided  for  both  crack  tips;  however,  in  this  problem  it  is 
assumed  that  the  crack  tip  at  the  rivet  hole  will  not  grow  through 
the  hole;  hence,  only  the  inside  tip  must  be  considered.  Determine  the 
crack  aspect  ratio  and  the  normalized  eccentricity: 


3.1875/2 


bx  1.5  + 1 + 0.1875/2 


= 0.61  = 


a + 1 + J 


e 

b 


12.0  - 2.59 


9.40 


= 0.78 


12.0  12.0 
From  Figure  7.39,  the  correction  factor  F is 
F = 1.11 

Let  critical  conditions  apply 


Kt  -»  K , a -»  a , a 
I c’  c 


Then, 


a = 
c 


a = 
c 


t) 

85  ksi  y/inT 
S/JT(  1.594  in.)  1.11 


a = 34.2  ksi 
c 
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Hence,  the  panel  is  fracture-safe  under  the  30  ksi  operating  stress 
as  long  as  the  other  end  of  the  crack  remains  stopped  at  the  rivet 
hole. 


The  easisst  way  to  approach  the  second  part  of  the  problem  is  to 
plot  the  stress-crack  length  history;  i.e.,  compute  and  plot  the  criti- 
cal stress  for  a number  of  assumed  crack  lengths.  The  results  are  shown 
in  Table  6.1  and  Figure  6.6. 

TABLE  6.1.  CRITICAL  STRESS  FOR  VARIOUS  CRACK  LENGTH  FOR 
PANEL  GEOMETRY  OF  FIGURE  6.5 


1.406 

1.906 

2.406 

2.906 

3.906 

4.906 

5.906 

6.906 

7.906 

8.906 


a 

e 

e = 

12  - bl 

e 

F 

m 

(in.) 

b 

0.36 

10.594 

m 

m 

0.52 

10.094 

m 

fm 

0.62 

9.594 

0.80 

1.12 

0.69 

9.094 

0.76 

1.16 

0.77 

8.094 

0.67 

1.24 

0.82 

7.094 

0.59 

1.32 

0.85 

6.094 

0.51 

1.39 

0.87 

5.094 

0.42 

1.46 

0.88 

4.094 

0.34 

1.52 

0.90 

3.094 

0.26 

1.59 

c 

(ksi) 


Figure  6.6  shows  that  up  to  a crack  length  of  about  4.0  inches 
the  panel  would  be  fracture  safe.  An  effective  crack  length  of  4 inches 
corresponds  to  a physical  crack  length  of  3.81  inches  growing  from  the 
edge  of  the  hole. 


* . ftr  . £ 


CRACK  LENGTH.  2a,  (in.) 

Figure  6.6.  Critical  fracture  stress  for  panel 
geometry  of  Figure  6.5. 

6.6  Crack  in  a Stiffened  Panel  [ 45] 

Figure  6.7  shows  a built-up  structure  consisting  of  a skin 
stiffened  by  a number  of  stringers  riveted  to  the  skin.  A crack  grows 
from  one  of  the  rivet  holes  as  shown.  The  fracture  stress  of  the 
panel  will  be  investigated  for  various  assumed  crack  lengths.  It  is 
assumed  that  the  9-inch  bars  are  repeated  and  that  the  crack  grows  in 
a symmetric  manner.  To  solve  this  problem,  it  is  necessary  to  separate 
it  into  two  parts: 

a)  A center-cracked,  edge-stiffened  sheet  with  the  center 
stiffener  removed. 

b)  The  effect  of  the  center  stiffener. 


Consider  the  center-cracked,  edge-stiffened  problem.  The  geometry 
of  Figir  e 6.8  applies.  In  the  solutions  of  Figures  7.65  through  7.69, 
the  edge  stiffeners  are  considered  as  integral ; however , it  is  felt 
that  the  rivets  introduce  negligible  error  for  crack  aspect  ratios  < 0.8. 
As  the  crack  approaches  the  edge  stiffeners,  the  influence  of  these 
stiffeners  predominates  and  the  crack  geometry  of  Figure  6.9  applies. 
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MATERIAL:  7075-T6 

WIDTH  BETWEEN  STRINGERS.  W:  9 in. 

RIVET  PITCH,  p:  1.00  in. 

SKIN  THICKNESS.  B:  0.063  in. 
STRINGER  THICKNESS:  0.260  in. 

(NO  CRACK  BUCKLING) 


Figure  6.7.  Stiffened  panel  with  a crack. 


110 


mt  * f*'  .< 


Figure  6.8.  Center -cracked  edge- 
stiffened  panel. 


Figure  6.9.  Crack  tip  approaching 
a stiffener. 

Ill 
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Step  1:  Determine  a,  the  nondimensional  inertia  parameter.  In 


this  case,  1^,  the  moment  of  inertia  of  the  edge  stiffener  (2-inch 


angle,  1/4 -inch  thick)  about  the  stringer  neutral  axis  perpendicular 
to  the  plane  of  the  skin,  is 


I = 0.347  in. 
z 


_ str 
3 E 


0.347 


0.063  (y) 


3 io7 


~ 0.01 


and  the  correction  factors  of  Figure  7.66  apply. 


Step  2:  Compute  the  extensional  stiffness,  S,  for  the  finite  width, 

cracked  strip.  From  Figure  7.66, 


A E 
str  str 


S-  0.25  (2  +.U7JD .IS.  = ,.65 


0.063  (-y)  107 


Step  3:  Prepare  a table  of  correction  factors  (\')  for  various 

crack  aspect  ratios  from  Figure  7.66  (Table  6.2). 


TABLE  6.2.  CORRECTION  FACTORS  FOR  THE  CENTER-CRACK, 
EDGE-STIFFENED  PANEL 


Total  Crack  Length 
2a  (in.) 


Crack  Aspect 
Ratio 
2a/w 


(From  Figure  7.66) 


h. it  .i 


i 


As  the  crack  approaches  the  edge  stiffeners,  2a/w>  0.8,  Figure  7.66 
does  not  apply.  The  problem,  however,  approximates  that  of  a non- 
symmetric  crack  approaching  a riveted  stiffener,  the  geometry  of 
Figure  6.9. 

Step  4:  Determine  the  correction  factors  as  the  crack  goes 

through  the  edge  stringer  areas  and  beyond  (Figures  7.63  and  7.64). 
Compute  the  stiffness  ratio  S'  and  the  ratio  of  crack  length  to  rivet 
pitch  for  crack  aspect  ratios  >0.8  (Table  6.3). 


TABLE  6.3.  CORRECTION  FACTORS  FOR  A CRACK  TIP 
PASSING  THROUGH  A STIFFENER  AREA 


2a 

w 

2a  B E 

b A E 

str  str 

(t) 

(Note  1) 

(Note  2) 

0.9 

1.09 

4.05 

0.7 

1.0 

1.21 

4.50 

0.45 

1.1 

1.33 

4.95 

0.28 

1.2 

1.45 

5.40 

0.32 

Note  1:  The  crack  is  assumed  to  grow  through  a rivet  hole, 

thus  effectively  doubling  the  rivet  pitch,  p'  = 2p. 

Note  2:  From  Figures  7.63  and  7.64. 

Figure  6.10  shows  the  correction  factors  for  various  values  of  2a/w. 

It  is  apparent  that  the  stiffening  effect  of  the  edge  stringers 
decreases  the  stress  intensity  factor  sharply  as  the  crack  tip  passes 
through  the  stiffener  area. 

Step  5:  The  second  part  of  this  problem  involves  accounting  for 

the  effect  of  the  center  stringer  on  the  stress  intensity  factor. 

This  can  be  done  by  using  the  solution  of  either  Figure  7.61  or 
Figure  7.62,  depending  on  whether  the  crack  passes  through  a rivet 
hole.  Figures  7.61  and  7.62  were  computed  for  a Poisson's  ratio  of 
0.33  and  a rivet-diameter-to-rivet-pitch  ratio  of  d/2p  = 0.1.  The 
correction  factor,  however,  is  insensitive  to  variations  in  this  ratio. 
Since  it  is  here  assumed  that  the  crack  grows  from  a rivet  hole, 

Figure  7.62  will  be  used.  The  resulting  correction  factors  are  shown 
in  Table  6.4. 

Step  6:  Multiply  the  correction  factors  in  Tables  6.2  and  6.3 

with  those  in  Table  6.4  to  obtain  combined  correction  factors,  \*  or  \**. 
These  are  shown  in  Table  6.5. 
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Figure  6.10.  Correction  factors  for  symmetrically  cracked 
plate  with  stiffened  edges. 


Step  7:  With  the  correction  factors  known,  it  is  now  possible 

to  solve  for  the  critical  stress.  One  must  have  a value  of  K for  the 

c 

particular  material  and  thickness  under  consideration,  For  illustration 

purposes,  a value  of  K =85  ksi  - Vin.  will  be  used.  The  critical 

c 

stress  is  computed  from 


a 

c 


K 

c 


(\*  or  \**) 


Results  are  shown  in  Table  6.6  and  in  graphical  form  on  Figure  6.11. 
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TABLE  6.6.  CRITICAL  STRESS  FOR  A 
STIFFENED  PANEL 


2a 

w 

a 

0 

c 

0.1 

0.9 

69.2 

0.2 

1.8 

52. 6 

0.3 

2.7 

42.9 

0.4 

3.6 

37.4 

0.5 

4.5 

33.7 

0.6 

5.4 

31.3 

0.7 

6.3 

29.4 

0.8 

7.2 

29.3 

0.9 

8.1 

35.8 

1.0 

9.0 

51.6 

1.1 

9.9 

80.2 

1.2 

10.8 

66.3 

It  is  seen  that  the  fracture  stress  decreases  until  the  crack 
approaches  the  edge  stringers.  At  this  point  the  edge  stringers  start 
to  pick  up  the  load  and  the  maximum  fracture  stress  reaches  the  material 
yield  stress  at  2a/w  ~ 1.05.  Hence,  the  stringers  should  be  effective 
in  arresting  the  growth  of  the  crack.  In  the  presence  of  the  panel 
stiffeners,  crack  buckling  should  be  minimal. 

6.7  Fatigue  Life  of  an  Edge  Crack 


An  edge  crack  of  1-inch  length  has  developed  in  the  20-inch-wide 
panel  as  shown  in  Figure  6.12.  The  plate  thickness  is  0.100  inch  and 
the  material  is  7075-T6  aluminum,  In  this  example  the  plane-stress 

1/2 

fracture  toughness  Kc  for  this  thickness  is  assumed  to  be  68  ksi-in. 

The  operating  cyclic  stress  of  the  part  is  sinusoidal  with  a maximum 
stress  of  25  ksi  and  a minimum  of  7.5  ksi  (R  of  0.3).  How  many  cycles 
will  be  allowed  before  critical  conditions  prevail? 

Write  the  stress  intensity  factor  as 

K = \a  yjfa 
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Figure  6.12.  Edge  crack  in  a panel. 


sl  j. 

Calculate  the  crack  aspect  ratio:  — = — = 0.05.  From  Bowie's  solu- 

tion (Table  3.4),  \ = 1.15.  Then 

K = 1.15  (25  ksi)  Vn(l)  in.  = 50.95  ksi-in.1^ 
o 

According  to  Forman  [66],  for  this  material  the  constants  n and  C are: 
n = 3 

-13 

C = 5 X 10 

From  Equation  (6.6), 

2 

N — N — . i o 2 2 

C ° (1.15)  «(5  X 10"  ) (1  - 0.3)  (25,000) 

68,000  ln  68,000 

50,950  1 50,950 

N - N = 144  cycles, 
c o 
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Tlius,  for  a sinusoidal  loading  of  25  ksi  maximum,  the  1.0-inch  crack, 
will  approach  critical  conditions  in  approximately  144  cycles.  This 
is  a rather  low  remaining  life.  A low  value  is  not  unexpected  since 

Ko  = 50.95  ksi  - / in.  approaches  the  critical  value  of  68.0  ksi  ,/in. 

The  above  is  only  anestimate,  since  the  geometric  correction  fac- 
tor \ changes  slightly  as  the  crack  grows,  and  the  above  applies  only 
if  \ is  a constant. 

Estimate  the  change  in  \ during  the  fatigue  crack  growth  by 

estimating  \ when  the  crack  reaches  the  critical  length  a . 

c 

K = f3a  Jna 
c v c 


or 


a 

c 


.2  n2 

\ o 


Jt 


a 

c 


( 1.15) 2 ( 25,000) 


1.78  in. 


a 178 

Crack  aspect  ratio  — = = 0.089.  From  Bowie's  solution  (Table  3.4) 

\ = 1.19.  Therefore,  changes  only  slightly  (from  1.15  to  1.19)  during 
the  crack  growth  from  a^  to  ac . In  light  of  material  variability,  the 

noted  results  of  144  cycles  are  probably  sufficiently  accurate  for  an 
engineering  estimate.  For  estimating  the  fatigue  life  of  flaws  in 
specimens  of  more  complex  geometry , one  should  use  a computer-assisted 
numerical  integration  as  in  the  report  by  Engle  [72], 


r*' 


Chapter  7.  COMPENDIUM  OF  STRESS  INTENSITY  FACTORS 


The  following  is  a collection  of  stress  intensity  factor  solutions. 
It  is  not  a comprehensive  collection.  In  the  interest  of  brevity,  solu- 
tions are  limited  only  to  those  judged  to  be  of  potential  usefulness  to 
the  designer.  Thus,  many  solutions,  though  elegant  and  ingenuous, 
thought  primarily  to  be  of  academic  or  research  interest,  are  not 
included  here. 

Solutions  included  are  from  several  sources;  however,  most  of  them 
have  been  included  in  other  collections.  Sources  of  stress  intensity 
factors  which  may  be  of  additional  usefulness  to  the  designer  are  ASTT1 
STP  381  (11],  the  handbook  by  Tada  [15],  the  book  on  three-dimensional 
solutions  by  Parmeter  [75],  and  the  design  manual  by  Wilhelm  [45].  Most 
of  the  solutions  included  here  have  been  included  in  one  of  the  above. 
Therefore,  most  of  the  references  given  in  the  following  will  refer  to 
one  of  the  above  source  books  rather  than  the  original  paper  which  con- 
tained the  solution. 

Because  the  solutions  included  here  were  taken  from  several  sources, 
the  notation  may  occasionally  vary  from  that  of  the  other  chapters. 

Where  departures  occur,  an  attempt  has  been  made  to  include  sketches  and 
explanations  with  the  solutions  explaining  the  meanings  of  the  various 
symbols . 

It  should  be  noted  that  in  two-dimensional  solutions  when  a concen- 
trated force  appears  in  the  solution,  the  force  should  be  taken  to  be 
force  per  unit  thickness,  i.e.,  lb/in.  Similarly,  a concentrated 
in-plane  moment  means  moment  per  unit  thickness,  i.e.,  inch- lb/inch. 


a tin  7 — r cot  7 
a cot  7 + r tin  7 


Figure 

with 


Tt 


tin  7 ■ \fr rs" 


ACCURACY:  EXACT 
REFERENCE.  [15] 

7.2.  A straight  crack  in  an  infinite  plate 
skew  loading. 
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ACCURACY:  EXACT 
REFERENCE.  [151 

Figure  7.6.  A straight  crack  in  an  infinite 
plate  with  concentrated  force  and  concen- 
trated moment  applied  to  one  surface. 


ACCURACY:  EXACT 
REFERENCE:  [15] 

Figure  7.13.  Two  semi-infinite  cracks  in  an 
infinite  plate  with  remote  moments. 
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ACCURACY:  EXACT 
REFERENCE:  [15] 

Figure  7.14.  Two  semi-infinite  cracks  in  an 
infinite  plate  with  skew  loading. 

\ 

*. 

1 

; 

1 

i 


i 


ACCURACY:  ±1%  [15] 

REFERENCE:  [15] 


Figure  7.16.  A single  crack  growing  out  of  a hole 
in  an  infinite  plate. 
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F (a/R)  ONE  CRACK 


F (a/R)  TWO  CRACKS 


a/R 

(UNIAXIAL 

STRESS) 

(BIAXIAL 

STRESS) 

(UNIAXIAL 

STRESS) 

(BIAXIAL 

STRESS) 

0.00 

3.39 

2.26 

3.39 

2.26 

0.10 

2.73 

1.98 

2.73 

1.98 

0.20 

2.30 

1.82 

2.41 

1.83 

a 30 

2.04 

1.67 

2.15 

1.70 

0.40 

1.86 

1.58 

1.96 

1.61 

o.&o 

1.73 

1.49 

1.83 

1.57 

0.60 

1.64 

1.42 

1.71 

1.62 

0.80 

1.47 

1.32 

1.58 

1.43 

1.0 

1.37 

1.22 

1.45 

1.38 

1.5 

1.18 

1.06 

1.29 

1.26 

2.0 

1.06 

1.01 

1.21 

1.20 

3.0 

0.94 

(X93 

1.14 

1.13 

5.0 

0.81 

0.81 

1.07 

1.06 

10.0 

0.75 

tt75 

1.03 

1.03 

OO 

0.707 

a 707 

1.00 

1.00 

REFERENCE:  [111 

Figure  7.17.  Either  one  or  two  cracks  growing  out 
of  a hole  in  an  infinite  plate. 
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ACCURACY:  3%  [15] 

REFERENCE:  [15] 

Figure  7.21.  A single  crack  growing  from  an  elliptical 
hole  in  an  infinite  plate. 


7 .2  Two-Dimensional  Solutions  for  Cracks  in  Semiinfinite  Plates 
(Figures  7.25  through  7.31) 
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ACCURACY:  K,,  K,,  WITHIN  ONE  UNIT  OF  LAST  DIGIT 
K| 1 1 EXACT  [15] 

REFERENCES:  [11.  15] 

Figure  7.25.  Edge  crack  in  a semi-infinite  plate. 
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ACCURACY:  K,.  K,,  BETTER  THAN  0.5% 
K|||  EXACT  [15] 
REFERENCE:  [15] 


Figure  7.26.  Edge  crack  in  a semi-infinite 
plate  with  concentrated  forces  at  the 
crack  mouth. 


146 


: V:  ’Vt'T  - » 


< .* 


0.2  0.4  0.6  0.8  1.0 

REFERENCE  1151 


Figure  7.34.  A center  cracked  strip  loaded 
with  concentrated  forces  on  the  crack 
surface . 
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K|  = o v'na  F (a/b ) 


□ — f 


AT  THE  CRACK  LOCATION 

2.  THE  EFFECT  OF  h/b  IS  PR ACTICALLY  NEGLIGIBLE  FOR  h/b  ?10 


REFERENCES:  [11,15] 

Figure  7.40.  An  edge-cracked  strip  loaded  in  tension. 
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F|U  <a/b)  = 

ACCURACY:  F(|  BETTER  THAN  1%  FOR  ANY  a/b 
F 1 1 1 EXACT  [15] 

REFERENCE:  i 15] 


Figure  7.43.  An  edge-cracked  strip  with  shearing 
and  tearing  forces  at  the  crack  mouth. 
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F (a/b,  h/bl 


» a/b 


ACCURACY:  BETTER  THAN  1%  [15] 
REFERENCE:  [15] 


Figure  7.48.  The  center-cracked  tension  specimen 
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CRACK  LENGTH  b/(R„  - Rd 
REFERENCE:  [75]  ° ' 


ntrcnunuL.  i / %jj 

| Figure  7.51.  A cylindrical  pressure  vessel  with  an  internal 

I longitudinal  crack. 
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ACCURACY:  BETTER  THAN  1%  [15] 

REFERENCE:  [15] 

Figure  7.56.  An  edge-cracked  disk  loaded  with  all-around  tension 
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ACCURACY:  EXPECTED  TO  BE  BETTER  THAN  2%  1151 
REFERENCE:  (15J 

Figure  7.57.  A center-cracked  disk  loaded  with  two  crack-ope 
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ACCURACY 
RE FERENC  I 
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K.  = u \/na  \ 


a - HALF  CRACK  LENGTH 
p=  RIVET  PITCH 
B = PLATE  THICKNESS 
E = MODULUS  OF  PLATE 
ESTR  = MODULUS  OF  STRINGER 
STRINGER  £ a$tr  = AREA  OF  STRINGER 


a/p  = 0.25 


< 0.6 


cr 

8 0.2  I 


2 (aBI  E 
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CRACK  ASPECT  RATIO,— 
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STRINGER  <£ 
a/p  = Q-25  T 
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CRACK  ASPECT  RATIO. — 
W 

REFERENCE:  [45] 


Figure  7.63.  Crack  near  a riveted  stringer 
(S'  = 0.5  and  1.0). 
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CORRECTION  FACTOR.  X 


a = HALF  CRACK  LENGTH 
p=  RIVET  PITCH 
B=  PLATE  THICKNESS 
E = MODULUS  OF  PLATE 
E$tr  = MODULUS  OF  STRINGER 


W 

REFERENCE:  (45] 


Figure  7.64. 


Crack  near  a riveted  stringer  (S' 


2.0). 
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CORRECTION  FACTOR.  X 


S 


CRACK  ASPECT  RATIO, 

W 

REFERENCE.  (451 

Figure  7.65.  Edge-stiffened  center-cracked  plate  (Q 
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WHERE.  B-  PANEL  THICKNESS 
lz  - MOMENT  OF  INERTIA  OF 

STIFFENER  ABOUT  NEUTRAL  AXIS 
PERPENDICULAR  TO  PANEL  PLANE. 


CRACK  ASPECT  RATIO. 


Figure  7.69.  Edge-stiffened,  center-i  acked  plate  (a  = 0.09) 


EXTENSION AL  STIFFNESS.  S 


Three-Dimensional  Solutions  for  Surface  Flaws,  Embedded  Flaws 
and  Through  Flaws  (Figures  7.71  and  7.91) 


ACCURACY:  EXACT 
REFERENCES:  [11,15] 


Figure  7.71.  Penny-shaped  crack 
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'A  E (k) 

-rr/2 
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1/4 


1 sin*  0 + — cos*  9 

I .2 


\/l  - k2sin2^  ds? 


k2  * 1 - b2/a2 


ki,max=  ki  w_  in/ 2,1 

2,1  /— 

K(  (a  = b)  = — \/n$ 

V 

Ki  (a  -*  <*)  * o y/nb 


a \/nb 

E (k) 


NOTE:  SEE  TABLE  7.1  FOR  VALUES  OF  E (k) 

ACCURACY:  EXACT 
REFERENCES:  115,  45] 


Figure  7.72.  Elliptical  crack  loaded  in  tension. 
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TABLE  7.1.  ELLIPTIC  INTEGRALS,  E(k)  AND  K(k) 


0.8192 

0.8660 

0.9063 

0.9397 

0.9659 

0.9848 

0.9903 

0.9945 

0.9976 

0.9994 

1.0000 


Not*:  The  definition  of  integral*  In  the  •v><rve 

tabulation  are: 


K(k)  - f -~4  , the  cos*  I 

/l  • k^  ain^o 

elliptic  integral  of  the  first  kind  and 


'l  * fc  sin  gdg  f the  cMSplet* 


elliptic  Integral  of  the  seconi  kind. 

If  J 


I MhQ 


(ft  *■  * 


elliptical  crack  in  a plate  subjected  to  uniaxial  tension  (b/a 


06  L 


Lptical  crack  in  a plate  subjected  to  uniaxial  tension  (b/a 


LEGEND 


REFERENCE:  [75J 


00 'Z 


Figure  7.84.  An  elliptical  crack  In  a plate  subjected  to  uniaxial  tension  (b/a  = 0.10). 
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NOTE  SOLUTION  IS  THE  SAME  AS  FOR  THE  CASE  OF 
ACRACK  APPEARING  SUDDENLY  I DUE  TO 
PROJECTILE  PENETRATION.  SAY ) IN  A STRESSED 
PLATE. 

REFERENCE:  (831 

Figure  7.9A.  Transient  stress  Intensity  factor 
(or  sudden  impact  load  applied  to  crack 
surface. 
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I’ STATIC  '"'I'STATIC 


•TIME  NORMAL  IZED  WITH  RESPECT  TO  TIME  REQUIRED  FOR  A 
DILATATIONAL  WAVE  TO  TRAVEL  ONE  CRACK  LENGTH 


REFERENCE  (84) 

Figure  7.95.  Transient  SIF  due  to  oblique  impact 
it  1 dilatationnl  stress  wave. 
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LIST  OF  SYMBOLS 


C,  cv  C’ 


e , c , e , etc 

xx  yy  xy 


f(§),  F(f) 


F - P - iQ 


8ir  8i ii 


“ 2(1  + •> 


k - b2/a2 


Crack  length;  half  crack  length 
Critical  crack  length 

Effective  crack  length 

Original  crack  length 

Area  of  stringer  cross  section 

Crack  spacing;  dimension  on  crack  surface; 
half  plate  width 

Plate  thickness 

Semi-major  axis  of  elliptic  flaw 
Constants 

Rivet  hole  diameter 
Eccentricity 
Strain  components 

Young's  modulus 

Young's  modulus  for  a stringer 

Complete  elliptic  integral  of  the  second  kind 

Correction  factors  for  stress  intensity 
factor  solution 

Concentrated  force,  Eq.  (3.21) 

Strain  energy  release  rate 
Critical  strain  energy  release  rate 

Strain  energy  release  rate  for  Modes  I,  11, 
and  111,  respectively 

Shear  modulus  of  elasticity 

Beam  depth.  Fig  3.10 

Moment  of  inertia  of  a stiffener  about  the 
stringer  neutral  axis  perpendicular  to  the 
plane  of  the  panel 

Parameter  in  elliptic  integral 


K 

K* 

K 


*i*  Ku*  Sn 

KIC 

Kiscc 


M 

n 

X 

X. 


P 

P 


*1 

9 


R 

R » 
R_ 


itin  *wr 


Stress  intensity  factor 

Obsolete  definition  of  stress  intensity 
factor,  K*  • K lyfrT 

Critical  stress  intensity  factor 

Final  stress  intensity  factor  in  fatigue 
analysis 

Original  stress  intensity  factor  in  fatigue 
analysis 

Conditioned  stress  intensity  (actor 

Stress  intensity  (actors  (or  mode*  1,  II, 
and  111,  respectively 

i'Unc  strain  critical  stress  Intensity 
(actor 

Threshold  stress  intensity  (actor  (or  stress 
cor risen  cracking 

Length  of  center  cracked  anJ  double  edge 
cracked  spec  men* 

Nut 

constant  In  fatigue  equations 
Xwmfcer  of  cycles  in  fatigue  analysis 
dumber  of  cycles  at  fatigue 

Original  number  of  cycles 

Rivet  pitch 

Concentrated  force;  percent  square  fracture 
Load  used  to  calculated  ^ 

Shear  stress  load,  fig.  l.lt 

Surface  flaw  shape  parameter , 

|t(»OlJ  - 0.112  U/lfV 

' ' 

Polar  coordinate*  centered  a*  the  c reck  tip 
Plastic  rone  sire;  hole  radius;  load  ratio 

Lord  ratio, 

Diameter  of  plastic  rone  for  "Wr  111 
solut ion 
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